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SomMAIRE. — Aprés avoir rappelé Véquation régissant la pénétration du rayonnement dans un milieu diffusant ef résumé la 


méthode de résolution approchée de Chandrasekhar, nous donnons des exemples d’application de cette méthode & une couche de 
brouillard et ala mer, et nous discutons les approximations faites. 


ZUSAMMENFASSUNG. — Die Gleichung fiir das Eindringen einer Strahlung in ein streuendes Medium wird aufgestellt und. 
kurz die Methode von Chandrasekhar fiir eine Naéherungslisung angegeben. Dann folgen Anwendungsbeispiele dieses Verfah- 
rens auf eine Nebelschicht und auf das Meer und schliesslich eine Diskussion tiber die erzielten Ndherungen. . 


SumMMARY. — The equation governing the penetration of radiation in a scattering medium is given and Chandrasekhar’s approxi- 
mate method of resolution is summarised ; this method is applied to a layer of haze and to the sea and the approximations are 


discussed, 


La pénétration et la répartition du rayonnement 
électromagnétique et tout spécialement des ondes 
lumineuses dans les milieux diffusants naturels 
(brouillard et brume, eau de mer) présente un grand 
intérét pratique ; étude théorique en est simple et 
aboutit a des équations bien connues, mais dont la 
résolution exacte a été jusqu’ici impossible. Le tres 
important travail de S. CHANDRASEKHAR [1] sur les 
problemes de transfert de rayonnement a apporté 
récemment de nombreuses méthodes nouvelles de 
résolution, exacte ou approchée, dont nous avons 
appliqué une partie a cette étude de la diffusion dans 
le brouillard et dans la mer [ 2}. 


I. Méthode. 1. Equation de transfert du rayonnement. 
— Nous considérons une couche plane et homogene d’un 
milieu diffusant caractérisé par le coefficient népérien 
d’absorption totale K, le coefficient d’absorption appa- 
rente par diffusion k = w) K oli w < 1, et le coeffi- 
cient d’absorption propre 8 = (1 — mo) AK. Une direc- 
tion sera repérée par langle 6 avec la verticale ascen- 
dante et l’azimut » ; nous posons p» = cos 0; soit H 
la profondeur comptée a partir de la surface supé- 
rieure de la couche recevant le rayonnement incident 
et + = KH la profondeur optique correspondante ; 
l’épaisseur optique totale de la couche est +, = K de. 
Nous écrirons le coefficient de diffusion dans la direction 
(u, @) faisant angle © avec la direction d’incidence 


IK 
(u', 9’) sousla forme 7 p(cos ©) = 7 Plu, oH.) (*); 


wT 
: K 

p (cos @) est la fonction de phase et la courbe ee (cos) 

en fonction de © est Vindicatrice de diffusion. 

Le soleil produit a la surface de la couche diffusante 
un éclairement = / sur un plan perpendiculaire a sa 
direction de propagation (uy, 9’); le ciel au-dessus de 
la couche a une luminance uniforme J). Soit I (7;.u, 9,) 
la luminance dans la couche, a la profondeur tz, du 


OPTICA ACTA 


rayonnement diffusé se propageant dans la direction 
(4, 9). 

L’équation de transfert a laquelle satisfait cette 
luminance s’établit facilement : 


dl ents 1 ANE po I 
(1) — ( f 9) = Lie a) ?) >. Lea ple, 3 Ua ?,) Fett w 


4 TE 


{ iy 7°27 ; f 
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Les conditions aux limites sont : 


(2) (10 ;4<0) =I, 
( L(t 5 4 > 0) =0(**), 


la deuxiéme condition se réduit a £ borné quand 7, 
devient infiniment grand. 


2. Résolution de Véquation de transfert. — L’indica- 
trice de diffusion doit étre représentée sous une forme 
mathématique et nous choisirons pour cela un déve- 
loppement en polyndmes de LEGENDRE : 

N 
Y t P C © A 
>», wPi(cos ®) ; 


t= 


(3) p(cos 9) = 


d’apres la propriété d’additivité de ces polynomes : 


(3°) P(e, e352 9) = 

Mula 

>, >; (2— 80m) ow)" P7" (v.) P7"(u') cos m (e' — ¢) 
l=0 m=0 


(*) On a évidemment la relation : 
A ES ; 
27 -— p(cos ©) sin @ dO =k =a K. 
4 
e/ 9 
(**) Cette condition correspond au cas simple d’un sol noir ; 


on peut la remplacer par des expressions plus compliquées 
pour tenir compte d’une réflexion du sol. 


9 J 


ou San = 1 Sim =:0, 
Som = O's mm se 0; 


ee ee) 
eta se (1+ m) ! 


des polynémes de LEGENDRE. 
Si on développe 
p 


et P;" sont les fonctions associées 


N 
> I" (<3 w) cos m (9'5 — 9), 
m=0 


Tw cre i) 


on peut séparer l’équation de transfert en N équations 


intégro-différentielles en /’’(z; ») ne contenant plus 
que les variables + et wu. 

Pour résoudre chacune de ces équations, CHAN- 
DRASEKHAR emploie la méthode des ordonnées dis- 
crétes, en remplagant l’intégrale par une somme d’un 
nombre fini de 2n termes suivant la formule de quadra- 


ais +n 
ture de Gauss f. 7th) du ay 


=), 


aj f(vj) ob wy 

J)=—n 

sont les zéros de Pon(u) ; on est alors ramené a résou- 

dre Nsystemes de 2n équations différentielles linéaires, 

ce qui ne présente aucune difficulté mathématique. 
La solution peut se mettre sous la forme : 


(9) I” (734) = 


+n ym eck N 
a 
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A=m “. \ 


k” sont | ines de I’é sristi 
“ es racines de l equation caractéristique (*) 


(*) Dans le cas m = 1 et m = 0, il y ala solution double 
k, = 0 et le terme correspond dans l’impression de Oe Ee) 


Aye t Altes ve 0 a 
se réduit al? + 1_9 [(1-—#)=+ ui] 


- LENOBLE 


[Oer. Acra 


et 1!” des constantes d’intégration que lon détermine 
par les conditions aux limites. 


3. Remarque sur Vétat de polarisation du rayonne- 
ment. —- Dans ce qui précéde, nous nous sommes inté- 
ressés au rayonnement global sans considérer son 
état de polarisation, or la diffusion modifie toujours 
plus ou moins la polarisation du rayonnement incident 
et il peut dans certains cas étre utile d’étudier théori- 
quement ce probleme [3]. Pour cela CHANDRASEKHAR 
représente une grandeur caractérisant Vénergie du 
rayonnement, par exemple la luminance, par une 
matrice de la forme : 


rf 
8) Ps tn ees dye 


V 


correspondant aux vibrations suivant 2 axes perpen- 
diculaires arbitraires / et r du plan d’onde ; ce rayon- 
nement de luminance totale J = J, + J, peut étre 
considéré comme la somme d’un rayonnement pola- 
risé elliptiquement de luminance 


I)? 2 U2 4 y2 eis 


, ou J, et Z, sont les luminances 


Thy es | ir 


caractérisé par langle ~ des axes principaux avec | 
et r et par l’ellipticité tg 8, et d'un rayonnement na- 
turel de luminance ['") = 7 — J!) ; on a les relations 


i,—I1, = I'?) cos 28 cos 2x , 
(9) U = I'?) cos 28 sin 2y , 
V =I") sin 2p . 
Le coefficient de diffusion s’écrira sous la forme 
K 
ie (cos ©) ot P (cos ©) est une matrice de phase a 
Tw 
4 lignes et 4 colonnes, traduisant l’action de la diffu- 
sion sur les différents termes de la matrice éclairement. 
En utilisant les matrices rapportées aux axes paral- 
leles et perpendiculaires au plan méridien (*), on 
obtient une équation exactement semblable a l’équa- 
tion de transfert habituelle ; la résolution en est évi- 
demment plus longue. 


II. Application a la brume. — La brume et le 
brouillard constituent généralement des couches 
d’épaisseur finie, ou Vabsorption propre est négli- 
geable pour le rayonnement visible et ultraviolet 
proche. 

Nous avons appliqué la méthode ci-dessus (approxi- 
mation d’ordre n = 3) A une couche de brume sans 
absorption propre (m, = 1), caractérisée par lindica- 
trice de diffusion (**) P (cos ®) = 1+ 1,73 P, (cos 0) 
+ P, (cos ©), d’épaisseur optique totale r,=4,606 (***), 


(*) Il est important de bien faire les changements d’axe 
necessaire pour exprimer P(u, @ 5 ues 

(**) Nous discuterons plus loin le choix de cette indicatrice. 

(***) Rappelons qu’il s’agit toujours de densités optiques 
neperiennes ; la densité décimale correspondante est 2. 
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Rayonnement 
vers le bas 


| vers le haut 


Fic. 1. — Répartition des luminances dans une couche de brume éclairée par un ciel uniforme 
(t, = 4,606 — p (cos ®) = 1 + 1,73 P, (cos + @) P, (cos ©) — sol noir). 
Les courbes représentent, en coordonnées polaires I(t ; %)/Zg en fonction de 0 dans un plan vertical quelconque ; 


1— t= 0 (Sommet de la couche); 2—t = 1,151 


au-dessus d’un sol noir. Les figures 1 et 2 montrent un 
exemple des résultats [4] ; on y a porté la lumi- 
nance J (a un facteur constant prés) en fonction de 
Vangle polaire 8, au sommet, au quart, au milieu, 
aux trois-quarts de la couche et au sol ; la figure 1 
donne //J, dans un plan azimutal quelconque pour la 
couche éclairée uniquement par un ciel de luminance 
I, ; la figure 2 donne //F’, dans le plan azimutal con- 
tenant le soleil et dans un plan perpendiculaire au 
précédent pour la couche éclairée uniquement par un 
soleil a4 60° du zénith produisant l’éclairement */ au 
sommet de la couche sur un plan perpendiculaire a sa 
direction. Dans le cas pratique, ov il y a simultanément 
les deux modes d’éclairage, le résultat s’obtient immé- 
diatement en ajoutant le terme dt au soleil et le terme 
du au ciel. 

Nous avons considéré dans ces calculs que le sol 
était noir ; si on admet qu’il présente une réflexion 
diffuse suivant la loi de Lampert avec lalbédo r, 
on peut ajouter un terme correctif Jp(+ ; vw, 9) a 
I(z;, 9) en évaluant approximativement la lumi- 
nance du sol 
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3 3— 7 = 2,303; 4— 7 — 3,454; 5 — + = 45606601). 


[L(r; 4, o) +I p(t3 4, ¢) uw duds ; 


n 


r 72m (70 
I =— 
Tw 
e 


De a 


Tp satisfait a une équation analogue a léquation de 
transfert avec les conditions 


(1 = 
Pes (ORS a) cae) 
[In (ase > 0) =4, 


Il ne semble d’ailleurs pas nécessaire de chercher 
un calcul tres précis de ce terme, tant que nous ne pos- 
séderons pas plus de données sur la réflexion des 
sols. 

Un calcul analogue au nétre a été fait par J. Oxr- 
VIER [ 5, 6] pour la méme couche éclairée par le soleil 
a 78° du zénith ; il serait souhaitable que ces calculs 
soient étendus a diverses couches et divers éclairages ; 
un bureau de calcul moyen pourrait certainement éta- 
blir un nombre de diagrammes important, couvrant. 
ensemble des cas usuels et fournissant des renseigne- 
ments intéressants pour la navigation aérienne. 
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vers le bas 


| Rayonnement 


a) 


i vers le haut 


0,1 
———S Es | 
Rayonnement 
vers le bas 
b) 
| vers le haut 
lia, 2. — Répartition des luminances dans une couche de brume éclairée par le soleil A 60° du zénith. 


(t, = 4,606 — p (cos ©) = 1 + 1,73 P, (cos @)+ P, (cos ®) — sol noir). 
Les courbes représentent en coordonnées polaires I(t ; u)/F en fonction de 0 


a) dans le plan vertical contenant le soleil ; 
b) dans un plan vertical perpendiculaire au précédent. 
1 — + = 0(sommet de la couche) ; 2— + = 1,151 ; 3 — +=2,3038 4 — + = 3,454;5—7 = 4,606 (sol). 
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400 450 500 550 Amt 600 
Fig. 3. — Répartition spectrale du rayonnement diffusé vers le haut parla mer 
(profondeur infinie. — 8 = @ eau distillée — k = 0,17 — ciel uniforme blanc) 


105 0 = 76/1. 4 Sa = 

III. Application a la mer. — L’eau de mer présente 
une absorption propre importante et trés variable dans 
le spectre ; mo est alors un paramétre caractéristique 
de chaque eau pour chaque longueur d’onde. Sauf au 
voisinage immédiat des cétes, ot le rdle du fond et 
Pétude de Véclairage qu’il recoit sont importants, 
on a généralement des profondeurs assez grandes pour 
considérer la mer comme un milieu diffusant de pro- 
fondeur optique infinie. 

Nous avons appliqué la méthode décrite précédem- 
ment a une mer de profondeur infinie, caractérisée 
par lindicatrice de diffusion 


p (cos @) = mo [1 + 1,73 Py (cos @) + Pz (cos ©) ] 


avec diverses valeurs de mp, et éclairée soit par un ciel 
uniforme, soit par le soleil 4 60° du zénith (dans lair). 

Si on considére une mer ayant Pabsorption propre de 
Veau distillée et une diffusion neutre telle que m)=0,65 
pour la longueur d’onde 500 mu, éclairée par un ciel 
blane de luminance J), la lumiére renvoyée vers le haut 
a la répartition spectrale que montre la figure 3 [7 | 
ou on a porté J(0 ; v)/Z, en fonction de la longueur 
d’onde \ pour un certain nombre de directions ; cette 
répartition, rapportée a l’étalon C de lumiere blanche, 
correspond a une longueur d’onde dominante de 
480 my environ, avec une pureté d’excitation de l’ordre 
de 65 % augmentant quand la direction d’observation 
se rapproche de la verticale ; cette mer apparaitrait 
bleue a l’observateur situé dans l’air, par suite de son 
minimum d’absorption propre dans le bleu ; bien 
entendu, l’existence d’une diffusion sélective peut dans 
la réalité modifier plus ou moins cette couleur. 

Les expressions de la luminance montrent que, aux 
grandes profondeurs, tous les termes de la somme 
d’exponentielles deviendront négligeables devant Ca, 
k, étant la plus petite solution de l’équation caracté- 
ristique pour m = 0 ; on atteint alors un état d’équi- 
libre ot la répartition spatiale de luminance ne varie 
plus avec la profondeur, est indépendante de la répar- 
tition de l’éclairage en surface, présente la symétrie 
de révolution autour de la verticale et ne dépend que 
de m, et de lindicatrice ; la figure 4 montre cette 
répartition d’équilibre en coordonnées polaires pour 
diverses valeurs de m, et avec l’indicatrice choisie ; la 


1(0 ;0 = 49°)/Ip. 


1(0 ; 8 = 21°)/Ip. 


| Rayonnement | 


vers le bas | 


[ vers le haut 


Fic. 4.— Répartition des luminances dans la mer al’équilibre. 
Indicatrice p (cos ©) = ap[1 + 1,73 P, (cos ®) + P, (cos O)]. 
Les courbes représentent, en coordonnées polaires, 
I(t ;)/I(e ; ©) en fonction de 6 dans un demi-plan vertical 
quelconque 


oe eka sop ye ale : S| 
sels tata pal sleet By = 0,66. 


6 oF 


figure 5 donne la variation de ky avec mo. La décrois- 
sance du rayonnement en profondeur est alors expo- 
nentielle avec le coefficient d’extinction ky K < K. 
L’expérience [8] semble confirmer ces résultats ; dans 
ce cas la mesure de l’extinction et de la répartition 
du rayonnement devrait permettre d’atteindre les 
coefficients d’absorption propre et de diffusion de l’eau 
de mer. 


IV. Discussion. 1. Précision de la méthode. — Les 
propriétés du milieu diffusant étant supposées con- 
nues, et son indicatrice représentée exactement par 
Vexpression employée pour p(v, 93 u’,¢?’), ce que nous 
discuterons plus loin, la seule approximation faite 
consiste dans l’emploi de la méthode des ordonnées 
discrétes, qui remplace une intégrale par une somme 
finie de 2 n termes. La discussion de l’erreur ainsi 
introduite est, dans sa généralité, un probleme mathé- 
mathique compliqué. Nous nous bornerons a donner, 
sur 2 exemples, une idée qualitative de la précision a 
attendre. Pour la diffusion primaire avec une indica- 
trice isotrope, il est facile de calculer exactement 

P41 
fo Pose) de 
ow 


et de comparer le résultat a 
+n 
Soa Ps 3); 
To 


des n = 3 lécart est inférieur a quelques %%, sauf 
quand 7 devient infiniment grand (+> 20), c’est-a-dire 
au moment ow la diffusion primaire est totalement 
négligeable. 

Aux grandes profondeurs, on a vu qu il s’établit un 
état d’équilibre régi par la plus petite solution k, de 
Péquation de transfert ; sur la figure 5 on a porté 
cette solution A, en fonction de w, pour Vindicatrice 
p (cos ©) =a, (1 + 1,73 P, (cos ©) + P, (cos 0) ] 
a divers ordres d’approximation n ; on peut dire que 
n = 2 donne des résultats corrects a Péquilibre pour 
Wo > 0,5 et n = 3 pour a > 0,3; plus la diffusion 
est faible par rapport a absorption propre, plus il 
faut la calculer dans un ordre d’approximation élevée ; 
mais il ne faut pas oublier ici que la limite w) = 0 ne 
présente aucun intérét pratique pour notre étude, 
puisqu’elle correspond au milieu non diffusant. 


2. Choix de l’indicatrice de diffusion. — Nous avons 
employé dans les calculs qui précédent une indica- 
trice représentée par une série de polynémes de 
LEGENDRE ; il nest pas exclu, bien entendu, que des 
méthodes analogues de résolution puissent étre appli- 
quées avec d’autres formes d’indicatrices. Néanmoins 
la propriété d’additivité des polynémes de LEGENDRE 
présente de tels avantages pour la simplification des 
calculs qu’il semble souhaitable de remplacer V’indi- 
catrice réelle du milieu par un développement en série 
de ces polyndmes. 


Indicatrices réelles. — \/indicatrice d’un milieu dif- 
fusant, caractéristique importante et physiquement 
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0,5} 
oS — 
OS @, 4 
Fig. 5. — Variation de k, en fonction de ap, a divers ordres 


d’approximation, pour Vindicatrice 
p (cos ©) = am [1 + 1,73 P, (cos @) + P, (cos ©)]. 


bien définie du milieu, est en général tres mal connue, 
son étude expérimentale directe est délicate au labo- 
ratoire et elle le devient encore beaucoup plus quand 
on veut travailler sur le terrain dans les milieux natu- 
rels comme la mer ou le brouillard. 

Dans le cas de particules diffusantes sphériques, 
présentant ou non de l’absorption propre, la théorie 
de Mie donne le coefficient de diffusion et la forme de 
Vindicatrice pour chaque dimension de particules ; 
connaissant la granulométrie du milieu, la théorie per- 
met donc de définir ses propriétés diffusantes ; pour 
les tres grosses particules, le calcul des formules de 
Mre devient tres pénible (*) et les lois de Poptique 
géométrique (diffraction et réflexion-réfraction) peu- 
vent étre appliquées. Pour la brume et le brouillard, 
ou les particules diffusantes sont des gouttelettes 
deau, ces calculs semblent donner des résultats en 
trés bon accord avec Pétude expérimentale directe. 
Pour la mer, nous possédons peu de renseignements 
sur la nature des particules diffusantes, mais il est 
probable quwil en existe de forme plus ou moins com- 
plexe et non sphérique ; dans ce cas, l’étude théorique 
devient difficilement applicable et la seule solution 
reste l'étude directe des propriétés diffusantes du 
milieu. 

D’une fagon générale, il semble que dans les milieux 
naturels (mer et brume) les particules diffusantes sont 
d’assez grandes dimensions par rapport a la longueur 


Pe Fae ie : , Q2ar 

(*) Le caleul a été fait (9) jusqu’a «a = Sao = 40, soit pour 
des diamétres de l’ordre de quelques microns pour le rayonne- 
ment visible. 
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d’onde et donnent des indicatrices de diffusion A maxi- 
mum tres marqué vers l’avant. 


Développement des indicatrices en polynédmes de 
LecenpRE. — Le développement en série de poly- 
nomes de ces indicatrices a rapide variation est évi- 
demment assez difficile ; avec peu de termes, toute 
formule, qui conserve bien le maximum avant, repré- 
sente mal la diffusion dans les autres directions et 
vice-versa ; il faut donc choisir un compromis, ou, si 
lon veut une représentation correcte, accepter de 
garder dans le développement un grand nombre de 
termes. 

Cuu et Cuurcnixy [10], développant une idée de 
HarrEL, ont montré récemment dans un tres inté- 
ressant travail que les formules de Miz peuvent faci- 
lement étre transformées pour étre mises sous la forme 
dune somme de polynémes de LeGENDRE et ils ont 
entrepris les calculs numériques qui donneront les 
coefficients de ces développements ; le nombre de 
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termes a garder augmente évidemment avec la dimen- 
sion des particules et serait de ordre de 2 «. Pour les 
tres grosses particules, il semble que l’on peut assimiler 
a la lumiere transmise la lumiere diffractée qui est 
concentrée dans un angle tres faible et utiliser pour 
la réflexion-réfraction un développement de l’indica- 
trice de WIENER. 

Pour représenter en série de polyndmes de LEGEN- 


K wear K - 
DRE 7? (cos ©) = iq 70 1+ >; Caer Moe -) 
n=1 \ 

une indicatrice donnée par ses valeurs numériques (indi- 
catrice de WIENER ou indicatrice expérimentale d’un 
milieu), diverses méthodes sont possibles : la plus 
élémentaire consiste a choisir arbitrairement N valeurs 
de © et a déduire les «, du systeme a N équations 
ainsi obtenu ; cette méthode donne évidemment des 
résultats trés variables suivant le choix des points et 
semble peu satisfaisante. On démontre mathématique- 


Fic. 6. — Développements de l’indicatrice de WIENER. 


—____—_——. p (cos @) d’aprés WIENER. 
as D (cos @) = 1 + ae Sas as Pe (cos ©). 
—- sO) = 1 Dob cos, & 

Ba eh t} * 3,31 P, (cos @). 
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= 1 + 2,33 P, (cos @) + 
3,31 P, (cos ®) + 2,89 Ps (cos ©) + 
+ 2,47 P, (cos @) + 1,79 P; (cos @) + 


a1. 90) P, (cos ©). 
3 


eo Ses Pp (cos 0) 


8 “e 


ment qu’on obtient l'approximation en moyenne (*) 
en calculant les coefficients «, par la formule 


_ e+1 r 
Kw, — — | —_ ~ p( (cos) P,, (cos @) d (cos ©) ; 
4n 2n+1 


nous avons calculé ainsi les 7 premiers coefficients pour 
Vindicatrice de veer We == 2,00 ey le 
Ode de 2.4) oe —— 1st og — 1 0.93: 
Une autre méthode as semble physiquement la plus 
logique consiste a écrire que la proportion de flux 
diffusé dans un nombre N d’angles déterminés est 
conservée ; lA encore il existe une arbitraire dans le 
choix des angles utilisés ; en divisant l’intervalle 0-7 
de variation de © en intervalles égaux, on obtient 
dans les premiers ordres pour l’indicatrice de Sek 
des coefficients assez proches de ceux donnés ci-dessu 

au contraire le développement d’ordre 2 (%; = 173. : 
a, = 1) employé dans les paragraphes précédents, 
correspond approximativement a la conservation des 

flux 0°-90 flux 0°-45° 


flux 909-1809 flux 135°-180° 
montre a titre d’exemples quelques développements de 
Yindicatrice de WIENER. 

En conclusion, il semble done que si l’on connait 
bien Vindicatrice d’un milieu et que l’on désire en 
tenir compte avec une bonne précision, il faut employer 
un développement d’ordre élevé, qui entraine des cal- 
culs longs, mais trés possibles avec les machines 
modernes ; si au contraire on veut arriver aux résul- 
tats par des calculs rapides, il faut se contenter, 
comme nous l’avons fait dans les paragraphes précé- 
dents, d’un développement d’ordre faible, qui ne peut 
représenter que grossierement l’indicatrice du milieu. 


IS Ca = 


. La figure 6 


rapports 


3. Influence de Vindicatrice de diffusion. — Un 
probleme qui se pose alors est de savoir si la forme de 
Pindicatrice a une influence importante sur la péné- 
tration du rayonnement et par conséquent si l’on ris- 
que de commettre une erreur appréciable en utilisant 
une indicatrice assez différente de lVindicatrice réelle 
du milieu. L’étude générale de ce probléme a partir de 
Péquation de transfert semble complexe ; mais on 
peut se donner une idée des phénoménes en compa- 
rant des résultats obtenus dans le méme cas avec des 
formes d’indicatrice différentes. Les figures 7 et 8 
donnent des exemples de cette comparaison ; la figure 7 
se rapporte a l’état d’équilibre qui s’établit dans un 
milieu de profondeur infinie, tandis que la figure 8 
concerne une couche de brume sans absorption propre 
d’épaisseur totale +, = 4,606 au-dessus d’un sol noir 
et éclairée par un ciel uniforme. Dans les deux cas, 
lécart des résultats n’est pas énorme pour deux formes 
d’indicatrices aussi totalement différentes que 


(*) C’est le développement qui rend minimum 
BAN g: 
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Fic. 7. — Etat d’ équilibre pour les indicatrices : 
—— p (cos®)=a [1+ 1,73 P,(cos@) + P, (cos@)] . 
-----=--- Pp (cos®) =a 

L,= I(u = 0 
ae 1) et log ee) 


Les courbes représentent k,, log G2 iGo 


en fonction de ap. 


p(cos ©) = to 
et 


p (cos 9) =w [1+ 1,73 P, (cos ©) + P, (cos ©) ] ; 


pour la brume les points calculés par OLivieR, avec 
une méthode différente, ont été obtenus en attribuant 
au milieu une indicatrice assez voisine de celle de 
GUMPRECHT pour x = 30 ; cette indicatrice a une 
allure analogue a 


p (cos ©) = wo [1+ 1,73 P, (cos 9) + P, (cos ®)], 


mais s’en distingue essentiellement par le maximum 
avant beaucoup plus important ; les résultats pour ces 
2 indicatrices sont tres voisins. Bien entendu, ’hypo- 
thése dune diffusion isotrope est trop éloignée de la 
réalité pour donner des résultats acceptables, mais il 
semble que, en général, il n’est pas nécessaire d’em- 
ployer rigoureusement lindicatrice réelle du milieu 
pour obtenir mieux que la précision de 10 % qui est 
tres acceptable dans l’étude des phénomeénes naturels. 

Le rayonnement venant d’une direction détermi- 
née ayant subi une seule diffusion a exactement la 
repartition donnée par l’indicatrice, alors que le rayon- 
nement ayant subiun grand nombre de diffusions sous 
des angles divers, doit étre tres peu influencé par la 
forme de l’indic atrice : il serait donc intéressant d’envi- 
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Rayonnement 
vers le bas 


—_——— 


( 


[vers le haut 


Fia. 8. — Répartition des luminances dans une couche de 
brume (a) = 1—t, = 4,606 — ciel uniforme — sol noir) pour 
les indicatrices : 

p (cos @) = a [1 + 1,73 P, (cos) + P, (cos ©)], 

= — = —_ p (cos) —a,; 

-------- indicatrice de GUMPRECHT pour «% = 30 (calculs 
de J. OLIVIER). 


La figure représente, en coordonnées polaires, J en fonction 
de 8 dans un demi-plan vertical quelconque, pour le sommet 
de la couche, le milieu et le sol. En trait épais luminance 
du ciel Jp. 


sager la proportion dans le rayonnement total des 
divers ordres de diffusion ; dans les cas ov la diffusion 
rimaire est importante et ot la couche est éclairée 
par le soleil direct, il est souhaitable et facile de faire 
une correction tenant compte de la forme de lindica- 
ce pour ce terme que Yon peut calculer exactement : 


RAYONNEMENT DANS LES MILIEUX DIFFUSANTS 
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4. Role des divers ordres de diffusion. — Si nous 
séparons la luminance globale J en une somme de 
termes correspondant a la luminance du rayonnement 
ayant subi 0, 1, 2... n... diffusions : ‘ 


(41) I(t; 2,9) = 2P%(c3-2,0) + (xs 0,9) +2. 
Se CES the) et are 


on peut remplacer l’équation de transfert par le sys- 
teme d’équations : 


Z ee Ya 
p SABE! = 12; u,9) 
di*(75 4,4 
o He a LAG 852) eee 
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POUT aioe: 
avec les conditions aux limites : 


\ P#9(0;n< 0) =Oet (0; p< 0) =] 


13 
) (Gr t= 0) — ©. 


On peut résoudre ces équations par itération, a partir 
de /1,en séparant les termes en cos m(g — °,) et en rem- 
placant Vintégrale de — 1a + 1 suru’ par une somme 
d’un nombre fini de termes, comme pour |’équation 
de transfert. 


Exemple. — Nous avons fait ainsi le calcul des 
divers ordres de diffusion jusqu’a 6, ainsi que le calcul 
du rayonnement global, dans le cas simple de la diffu- 
sion isotrope (p(u, 9; uv’, 9’) =m) pour une couche de 
profondeur infinie, éclairée uniquement par le soleil a 
60° du zénith. 

La figure 9 montre les résultats obtenus sous la 

3 6 
forme /1/J, > T/T et » I"/Ten fonction de + et pour 


n=1 n=1 
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Fic: 9. — Role des divers ordres de diffusion. 
Couche de profondeur infinie — diffusion isotrope — soleil a 60° du zénith — en haut a) = 0,5 ; en bas a) = 1. 
Les courbes représenten t 


3 6 
PI, S T/T et » I4/I 


n=1 
en fonction de t pour les directions 
0 = 1599 ——___ 
) ee Oe ae 


une série de directions, dans le cas m)= 1 (milieu sans 
absorption propre) et dans le cas m) = 0,5 (absorp- 
tion propre égale a l’absorption par diffusion). Bien 
entendu, ces courbes doivent varier avec la forme de 
Pindicatrice choisie, les conditions d’éclairage, et étre 


n=1 
0 = 760 eee 
@= 210 


modifiée par la présence d’un sol a profondeur finie, 
mais elles nous semblent intéressantes par les conclu- 
sions qualitatives qu’elles permettent de tirer. La 
part des diffusions d’ordre faible et en particulier de 
la diffusion primaire est minimum pour le rayonne- 
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ment de bas en haut et maximum pour le rayonne- 
ment a la verticale de haut en bas, ce qui s’explique 
facilement par le réle du rayonnement diffusé pres de 
la surface et transmis par le plus court chemin a la pro- 
fondeur considérée. Dans les milieux sans absorption 
propre, les diffusions d’ordre élevé deviennent tres 
rapidement prépondérantes ; dans notre exemple la 
diffusion primaire ne représente que 10 °% du rayon- 
nement global dés une profondeur optique de l’ordre 
de + = 0,5 pour le rayonnement de bas en haut et de 
Pordre de + = 1,5 pour le rayonnement de haut en 
bas ; 4 t ={4, les 6 premiers ordres de diffusion repré- 
sentent pour toutes les directions moins de 10 % du 
rayonnement global. Quand il y a un terme d’absorp- 
tion propre, une méme profondeur optique corres- 
pond a un moins grand nombre de particules diffu- 
santes et Vimportance des diffusions d’ordre élevé 
naugmente que beaucoup plus lentement ; pour 
Bo — 0,5 toujours dans notre exemple, la diffusion 
primaire est réduite a 10 °%% du rayonnement global 
vers + = 3 pour le rayonnement de bas en haut et 
seulement vers + = 15 pour le rayonnement de 
haut en bas ; pour cette direction de haut en bas, les 
6 premiers ordres représentent pratiquement tout le 
rayonnement jusqu’a +t = 4, les 3 premiers ordres 
fusquart =4" 

Il semble donc que pour les milieux ot la diffusion 
est prépondérante, dés une faible profondeur, le cal- 
cul global avec une indicatrice de diffusion grossiere 
doit donner des résultats raisonnables ; au voisinage 
immédiat de la surface (*) et surtout pour le rayonne- 
ment dirigé vers le bas, il sera bon de calculer exacte- 
ment le terme de diffusion primaire avec l’indicatrice 
réelle et de faire la correction. Pour les milieux, oti 
Vabsorption propre est importante, au voisinage de la 
surface, il sera souvent préférable de calculer exacte- 
ment par itération les premiers ordres de diffusion 
qui représentent la presque totalité du rayonne- 
ment (**); quand la profondeur optique augmente, le 


(*) Par une discussion un peu différente, OLIVIER a mon- 
tré également qu’au voisinage de la surface, il faut tenir 
compte de la forme exacte de l’indicatrice et surtout de son 
important maximum avant. . 

(**) Généralement dans ce cas, le rayonnement décroit 
approximativement en progression géométrique quand l’ordre 
de diffusion croit et on peut utiliser cette loi pour évaluer lc 
petit terme correctif a ajouter aux premiers ordres calculés. 
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calcul global redevient le plus avantageux, mais il sera 
bon de calculer exactement le terme de diffusion pri- 
maire jusqu’a des profondeurs assez importantes. 


V. Conclusion. — La méthode de CHANDRASEKHAR 
permet par des moyens mathématiquement simples, de 
résoudre l’équation de transfert qui régit la pénétra- 
tion du rayonnement dans un milieu diffusant ; en 
employant des développements a un nombre suffisant 
de termes, ce qui est possible avec les machines & cal- 
culer modernes, il semble que Von peut atteindre 
toute la précision désirée. Si plus modestement on se 
limite aux calculs possibles a la main, avec des déve- 
loppements a 2 ou 3 termes, on obtient sans peine des 
résultats, qui dans l’ensemble sont certainement cor- 
rects a mieux de 10 % et donc parfaitement utili- 
sables dans la pratique (visibilité dans la brume, éclai- 
rage sous-marin) ; mais alors, comme toute approxi- 
mation, ces calculs doivent étre appliqués sans perdre 
de vue la réalité physique et il peut étre nécessaire 
dans certains cas de les compléter et de les corriger 
par une évaluation précise des premiers ordres de 
diffusion, ceci en particulier quand la couche est 
éclairée par le soleil, au voisinage immédiat du som- 
met sil n’existe pas d’absorption propre, et jusqu’a 
des profondeurs d’autant plus grandes que le terme 
d’absorption propre devient plus important par rap- 
port a la diffusion. 
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Image quality and optical resolution 


E. H. Lixroor, The Observatories, University of Cambridge 


SumMmARY. — A brief account is given of the reasons why the classical nolion of resolving power fails lo provide a satisfactory 
qualily criterion of optical systems designed to form images of extended objects. The discussion indicates that a quan- 
tity here called the structural resolving power, which is formally identical with an assessment proposed independently by 
SCHADE (1952) and FELLGETT (1953), seems well suited to replace the old notion of resolution limit in evaluating the tmages 
of extended objects. On the basis of known results about images by systems with annular apertures, and of the rather close 
relation between the structural resolving power of an optical system and its information-passing capacity, il ts concluded that 
the image of an extended object showing low-contrasl fine structure begins to show vistble deterioration when the reduction in 
the structural resolving power of the optical system by aberrations or by central obstruction reaches about 20 per cent. 


Sommaire. — L’auleur rappelle briévement les raisons pour lesquelles la nolion classique de pouvoir séparateur ne peut fournir un 
criterium de qualité satisfaisant pour les systémes optiques calculés pour la formation @’ images @ objets étendus. el} 

La discussion montre qu’une grandeur appelée ici pouvoir séparateur de structure, fout a fait analogue a une grandeur définie 
indépendamment par ScuHADE en 1952 et par FELLGETT en 1953, semble bien convenir au remplacement de Vancienne notion 
de limite de résolution pour Vappréciation des images d objets étendus. 

En s’appuyant sur les résultats connus concernant les images formées par des systémes a ouvertures annulaires et sur la relation 
étroite existant entre le pouvoir séparateur de structure d’un systéme oplique et sa capacité de transmission d’information, on 
conclut que Vimage d’objets étendus ayant une structure fine faiblement contrastée commence a s’altérer visiblement quand la 
réduction du pouvoir séparateur de structure du systéme due aux aberrations ou a une obturation centrale atteint environ 20 %. 


ZUSAMMENFASSUNG. — Es werden kurz die Griinde dargelegt, warum der klassische Begriff des Auflésungsvermdgens als Giitekri- 
terium versagt, wenn man ihn auf optische Systeme anwenden will, die fiir die Abbildung ausgedehnter Objekte bestimmt sind. 
Eine Grosse, die hier “ strukturauflésung ” genannt wird und im Wesentlichen mit der von SCHADE (1952) und FELLGETT 
(1953) vorgeschlagenen Festsetzung tibereinstimmt, scheint jedoch gut geeignet zu sein, fiir die Beurteilung der Abbildung aus- 
gedehnter Objekte die alte Definition des Aufldsungsvermégens zu ersetzen. Die bekannten Ergebnisse bei der Abbildung durch 
Systeme mit ringformiger Pupille und die enge Beziehung zwischen der Strukturaufldsung und der Informationskapazitat optis- 
cher Systeme lassen schliessen, dass das Bild eines ausgedehnten Objektes mit feinen kontrastarmen Einzelheiten eine deutliche 
Verschlechterung erst dann zeigt, wann die Strukturauflosung des optischen Systems infolge von Aberrationen oder von einer 


zentralen Abblendung um 20 % zuriickgeht. 


I. Generalities. — The classical notion of resolving 
power was never precisely formulated except in a few 
very special cases, and has consequently been found 
insufficient for the needs of present-day instrumental 
optics. Something like a crisis has been reached in the 
last few years and, as stated by G. Toratpo in a 
recent paper [ 1 |, the theory of optical resolving power 
is at present undergoing a transformation. The diffi- 
culties which have hitherto prevented the idea of opti- 
cal resolution from receiving a general analytical for- 
mulation seem to arise in part from its having been 
too closely linked with the notion of a sufficient sepa- 
ration between the central peaks of the diffraction 
images of two neighbouring bright points in the object 
surface. This may have led to the assumption that 
the power of an optical system to reproduce fine detail 
in the images of continuous objects could somehow be 
measured by a single number, called the resolution 
limit. There is, for good reason, no generally agreed 
definition of the meaning of this term. The resolution 
limit of an aberration-free system is often taken to be 
that angular separation of two bright object points 
which places the geometrical image of each on the 
first Airy dark ring in the diffraction image of the 
other, and the ‘“* RAY Lercu resolution limit ” defined 
in this way agrees moderately well with the value 
(about 0”.57 of arc for an 8-inch aperture) adopted by 
Dawes as a result of his observations of double stars. 
But, besides involving rather artificial assumptions 
about what is meant by fully resolved, the very con- 
cept of a resolution limit presupposes, as was pointed 
out by Tora.po (I. ¢. [1 ]), that the observer knows in 


advance a great deal about the structure of the object ; 
in the present instance that it consists of either one 
or two bright points. Toraldo has further pointed out 
that, since many different intensity distributions in 
the object surface may yield exactly the same diffrac- 
tion image, it is only when the observer is able to 
choose from among these distributions on the basis of 
prior information about the object that it is permis- 
sible to speak at all of resolution in the conventional 
sense. 

In the majority of practical cases, the available 
prior information about the object is quite insufficient 
to make it possible to proceed in this way, and the 


- question arises whether an acceptable definition of 


resolution can be formulated which is applicable to the 
imaging of continuous objects about which little is 
known beforehand. The discussion of this question is 
the main object of the present paper. 

In many present-day applications the leading requi- 
rement is that previously unspecified fine detail in a 
continuous object shall appear in its image with suffi- 
cient contrast to emerge from the “ noise ” of a pho- 
tographic negative or print and with sufficient fidelity 
to be interpretable. Metamorphosis of fine detail will 
not necessarily prevent its being correctly interpreted 
by a practised observer, but it will in general make 
recognition more difficult and so reduce the total 
amount of detail recognised. 

Recognition of metamorphosed details can be regar- 
ded as a decoding of the image, and if this decoding 
entailed no loss of information in practice it would be 
by the power to discriminate, in its images, between 
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one object and another that an optical system could 
most realistically be assessed. As already mentioned, 
decoding by image interpretation is probably not 
fully efficient, and the use of image-reconstruction 
methods (optical or electronic) to improve fidelity is a 
complication which is not acceptable in many appli- 
cations. Nevertheless, it seems a reasonable starting 
point to consider.how far the notion of optical resolu- 
tion may agree with that of the discriminating power 
of an optical system, here defined as the logarithm of 
the effective number of distinguishable images [ 2 ] 
which it forms when presented with a statistically 
characterised object set. 

Discriminating power in this sense agrees in value 
with the statistical mean information content of an 
image for the given object set ; and its value depends 
on the characteristics of the object set as well as on the 
optical design. This dependence is entirely in accord 
with intuitive ideas about discriminating power and 
about information, but does not agree well with the 
conventional notion of optical resolution, which is 
usually thought of as a property of the optical system 
per se. It follows that we cannot obtain a satisfactory 
analytical realisation of the conventional notion of 
optical resolution simply by identifying it with discri- 
minating power. 

The essential difficulty can be traced to the fact that 
in order to specify the discriminating power of an opti- 
cal system for all object sets we require, not a single 
number, but a function, for example the 7 function 
(response function) of the system [3 ], while the notion 
of optical resolution purports to assess discriminating 
power, in some sense of the term, by a single figure of 
merit. 


II. Structural resolution. — A useful figure of merit 
can, however, be obtained by defining the optical reso- 
lution of a system as its normalised discriminating 
power for those object sets which, in a certain sense, 
test it most searchingly, namely random low-contrast 
object sets | o(x, y) | in which the structure is only just 
discernable through the photon noise [4]. It can be 
shown that the assumption of randomness (interpre- 
ted as greatest entropy/power ratio) implies that the 
statistical mean structure-power is equally distributed 
among all the spatial ‘“‘ sampling frequencies ” (Ws %q) 
which are passed by the system, that the different fre- 
quencies (w,,, %,) are statistically uncorrelated, and 
that the probability distribution of the coefficient «,,, 
of a given term «,,, exp | 2 wi (uw, x + 9, y) { in the 
double FourrerR expansion 


a(x, y) = > Xoq EXP 2 wi (u, Ps nae y)} (2; y) nA. 


Po 


‘atic 


of o(x,y) over the isoplanatism patch 


(x, y) notin A 


A(—astsa,—>b0e ys 5) 


is Gaussian to a sufficient approximation. 
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Then the discriminating power of the system wor- 
king in the isoplanatism patch A is given by the expres- 
sion [5 } 


Pein 
>, ee 2 ) 
(2iy eA 5 Te | t,(u, 9) |? du de, 


where 7, (u,v) is the response function (FourRIER trans- 
form of the “ spread function ”’) of the system in the 
isoplanatism patch A, while P and N denote respecti- 
vely the statistical mean structure power in the object 
set and the photon noise-power level, so that P/2N is 
proportional to the mean brightness. We normalise 
the discriminating power (2) by removing the factor 
{[A| P/2N and define the structural resolution or struc- 


tural resolving power of the optical system as the num- 
ber 


(3) lite = i| [z,(u, ¢) ? dude. 


It is very interesting that this quantity R,, rather 
closely related in the low-contrast case to the informa- 
tion-passing capacity of the system, is identical [ 6 | 
with an assessment of image quality already proposed 
independently by O. H. Scuape [7] and by P. B. 
Fre_icerr [8]. In a correctly focussed aberration- 
free monochromat with circular aperture, it is equal 
to 0.168 times the number of Arry discs in the recei- 
ving surface per unit solid angle in the object surface. 
The effect of defocussing is shown in figure 1, which 
was computed from H. H. Hopxins’s response func- 
tion curves [ 9 ]. It will be seen that even 1/4 fringe of 
defocussing reduces the value of R, by more than 
30 per cent. 


OS 1:0 fringes of 
defocussing 
Kia. 1. — Effect of defocussing on the structural 


resolving power of an aberration-free optical system. 
The dotted line corresponds to the Rayleigh quarter- 
wave limit. The ordinates show the values of RIK» 
where Ry = (/* sin? «)/0.696 22 is the structural resol- 
ving power of a perfectly focussed system. 
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Since there is some reason (see § 3 below) to believe 
that a barely visible deterioration in image quality 
corresponds to a drop of about 20 per cent in structu- 
ral resolution, figure 1 is in good agreement with 
Conrapy’s remark [10] that although the RAYLEIGH 
quarter-wave tolerance limit is abundently justified 
so far as the resolution of point sources is concerned, 
it might with advantage be tightened to 1/6 or even 
1/8 of a wavelength for instruments intended to form 
high-quality images of extended objects. 

In the practically important case of a non-isopla- 
natic system imaging polychromatic light, the local 
structural resolution at the point (x, y) of the field may 
be defined as the quantity 


(ay gees | 


where +(w, ¢; x, y) is the “ polychromatic local response 
function” at (a, y). 

The structural resolution AR of the system over a 
given working field / is then defined as the mean 
value of R,.,, over F : 


(5) R= van / | sp Ray dedy. 


The equivalent band radius r of the optical system 
working over the field # may be defined as the radius 
of a circle of area R in the frequency plane ; that is to 
say, by the equation 


wipe == Sate 


|r(u, 9; x, y)P du de, 


In an aberration-free monochromat the reciprocal 
diameter 1/27 of this circle is equal to 1.22 times the 
Ray Leicu resolution limit. The reciprocal diameter 
of the circle of spatial frequencies passed by the sys- 
tem is 7/7.66 = 0.41 times the RayLeran resolution 
limit. The radius of this circle defines the spatial fre- 
quency cut-off of the system, beyond which 7 is 
always zero; it corresponds to a line spacing in the 
image surface equal to 4 times the focal ratio. Aber- 
rations, including defocussing, leave the spatial fre- 
quency cut-off unaltered, but reduce the equivalent 
band radius by reducing | <(u, ¢) |. 

The content of these remarks may be summarised 
as follows : an optical system imaging incoherent 
objects is characterised by its monochromatic local 
response function +(w,v ; 2,y; 2). To abstract a figure 
of merit from this characterisation is equivalent, gene- 
rally speaking, to forming a weighted functional mean 
of +(w,0 ; vy ; A), and the present discussion indicates 
that the functional mean which seems to correspond 
best with the old idea of resolving power in polychro- 
matic light of given spectral distribution is the struc- 
tural resolution 


(6) f= all, avdy ff | <(w,9; x,y) 2 dude, 


where the polychromatic local response function 
t(U, 732, y) satisfies the equation 
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(7). xe Yar) = [ec ' L,Y) aU, 03H, Ys ry dr 


and (A; 2, y) is the effective spectral distribution of 
the light issuing from the object surface element 
situated at (2, y). The local structural resolution (4) is 
identical with an assessment of image quality already 
proposed by Scuape and Frxtieerr independently. 


Ill. Image quality tolerances in terms of structural 
resolution. — Some idea of the visual effect of a speci- 
fied drop in structural resolution may be gained by 
examining, as a particular case of practical importance, 
the reproduction of fine detail in paraxial images for- 
med by ordinary Newtonian reflecting telescopes. 
Astronomical observers agree that with an obstruction 
ratio of less than 0.25 there is no appreciable deterio- 
ration of the image, but that the use of a diagonal flat 
of (linear) obstruction ratio greater than 0.25 noticea- 
bly impairs the observation of planetary detail ; and 
similar conclusions have been arrived at by visual 
users of reflecting microscopes [ 12 }. 

The effects of central obstruction on the transmis- 
sion function of an aberration-free monochromat are 
easily calculated, since the response function 7(u, 9) 
is the area common to two annuli in the (uw, #)-plane 
(frequency plane), each similar to the clear aperture 


(see fig. 2), with their centres ata distancew = V u? + 2 


YS 


Fic. 2. — Calculation of response function for an 
aberration-free system with central obstruction. The 
value of t(u, v) is given, apart from a normalising 
constant, by the area of the shaded region when 
C, Cx Vu2 + v2, The radius of each of the large 


circles is a , Where 2 w, is the spatial cut-off fre- 
quency. 


apart. Figure 3 shows t-curves, obtained (with one 
exception) by measuring these areas with a planime- 
ter, for a selected set of values of the obstruction 
ratio [13]; the parameter S is defined as the ratio 
«/w», Where 2, is the cut-off frequency of the system 
for self-luminous objects. Figure 4 shows the effect of 
central obstruction on the structural resolution of the 
system. From this figure we see that a central obstruc- 
tion-ratio of 0.25 corresponds, in an aberration-free 
monochromat, to a drop of about 18 per cent in the 
structural resolution. 
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Fre. 3. — Response functions t(u, v) = t(s) of aberration-free monochromats with different amounts of 
central obstruction. The parameter s = @/a , Where wo = Vu -+- v® and 2 q@, is the spacial cut-off frequency 


for self-luminous objects. 
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central obstruction 
ratio 


The validity of the above conclusion is not restric- 
ted to monochromats, it holds without much change 
for aberration-free systems working in polychromatic 
light such as sunlight. In figure 5 are compared the 
response functions of an aberration-free system 
without central obstruction and of one with central 
obstruction 0.25, both in monochromatic light and 
in sunlight [14 ]. On calculating from equation (3) the 
value of #, corresponding to each of these four curves, 
it appears that a central obstruction of (linear) ratio 
().25 reduces the structural resolution of the system by 


OPTICA ACTA 


Fic. 4. — Effect of central obstruction on the 
structural resolving power of an aberration-free mo- 
nochromat. Ry = (f? sin? «)/(0.696 22) is the structu- 
ral resolving power of the unobstructed system. 


a factor 0.84 in sunlight, as compared with 0.82 in 
monochromatic light. 

As already mentioned, the structural resolution of 
an optical system is rather closely related to its infor- 
mation-passing capacity ; indeed the two become equi- 
valent when the system is used on a random low- 
contrast object set. Furthermore, a slight decrease in 
image quality shows itself mainly through its effect 
on the reproduction of low-contrast fine detail. It fol- 
lows that the conclusion drawn above, viz. that a 
decrease in R, of about 20 per cent corresponds to a 


16 E. i. 


0:2 0-4 0-6 0-8 


Fie. 5. 


LINFOOT 


— Response functions of an aberration-free system for “ grey ’ 
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object structure in monochro- 


matic light and in sunlight, with and without a central obstruction of linear ratio 0.25 ; (1) monochromatic 
light, no central obstruction ; (2) monochromatic light, central obstruction 0.25 ; (3) sunlight, no central 


obstruction; (4) sunlight, central obstruction 0.25. The parameter s = @/@) Where w = Vue + v2 and 2 w? is 


the spatial cut-off frequency at the peak wavelength Ag 


hardly perceptible deterioration in the image, may be 
expected on these general grounds to hold without 
much change for image-deterioration caused by small 
amounts of malfocussing, or of spherical aberration, 
or of aberrations of more general type, as well as for 
image-deterioration caused by central obstruction of 
an aberration-free aperture. 


Manuscrit reeu le 15 aott 1956. 
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Der Einfluss des Offnungsfehlers auf die Bildgiite 


Joacuim FockE 
Mathematisches Institut, Leipzig 


ZUSAMMENFASSUNG, — Die Lichtverteilung im Bild von Kante, Streifen und Gitter wird unter Beriicksichtigung des Offnungsfehlers 
auf wellenoptischer Grundlage theoretisch berechnel und sein Einfluss auf die Bildgiite an Beispielen gezeigt. 


SOMMAIRE. == IL) répartition de la lumiére dans Vimage de tests (bord, bandes ou grille) est calculée théoriquement en fonction de 
Vaberration sphérique, & partir de optique ondulatoire et son influence sur la qualité de Vimage est montrée dans plusieurs 


exemples. 


SUMMARY. — The distribution of light in the images of edges, strips and gralings in the presence of spherical aberration is cal- 
culated using wave theory, and its influence on image quality is illustrated by examples. 


Zur Beurteilung der Abbildungsgiite photographi- 
scher Objektive gibt es auch heute noch kein verbind- 
liches Kriterium. Die oft geforderte rein zahlenmi- 
ssige Leistungsangabe etwa in Form des Auflésungs- 
vermoégens hat sich als irrefiihrend erwiesen, denn 
abgesehen davon, dass sich dieses kaum absolut be- 
stimmen lasst, spielt bei der photographischen Bild- 
giite nicht das Auflésungsvermogen sondern die Kon- 
trastwiedergabe die ausschlaggebende Rolle. Man wird 
deshalb nach wie vor die Abbildung geeigneter Test- 
objekte wie Kante, Streifen, Kreisscheibe (1) unter- 
suchen miissen. Neuere Vorschlage gehen auch dahin, 
die Abbildung periodischer Objektstrukturen, vor 
allem Gitter mit sinusférmiger Helligkeitsverteilung, 
zu betrachten [1]. Da das Bild eines solchen Gitters 
stets wieder die gleiche sinusformige Intensitatsver- 
teilung aufweist, lasst sich die Abbildung durch nur 
einen Zahlenwert, die sog. Kontrastiibertragung (con- 
trast transmission), vollstiindig kennzeichnen, welche 
nur von der Linienfrequenz des Gitters abhingt. Man 
kommt somit obiger Forderung nach zahlenmassiger 
Leistungsangabe niher. Jedoch ist auch dieses Verfah- 
ren nicht frei von Fehlschliissen, da es teilweise 
«falschen Kontrast » anzeigt, welcher nur auf der 
periodischen Struktur des Objektes beruht und 
nicht allgemein auf die Bildgiite schliessen lisst. 
(Siehe hierzu [1], Figure 9 und Figure 10 und auch 
unsere Abbildung 3). 

In der vorliegenden Arbeit soll nun an Hand der 
genannten Teste rein theoretisch auf Grund wellen- 
optischer Rechnungen der Einfluss des Offnungs- 
fehlers auf die Bildgiite untersucht werden. Wir 
schliessen uns eng an die vor kurzem in dieser Zeit- 
schrift erschienene Arbeit des Verfassers [2] an, in 
welcher die Abbildung eines Objektpunktes bei vor- 
handenem Offnungsfehler eingehend theoretisch unter- 
sucht worden ist. 

Die Intensititsverteilung in der Punktkaustik 
kann somit als weitgehend bekannt angesehen werden. 


(2) Siehe hierzu den Vorschlag von J. HARTMANN auf der 
57. Tagung der Deutschen Gesellschaft fiir angewandte 
Optik, 1956. 


Hieraus ergibt sich dann die Lichtverteilung im Bild 
eines inkohirent strahlenden Objekts durch entspre- 
chende Uberlagerung. Speziell fiir die leachtende Halb- 
ebene (Kante), ihr geometrisches Bild in der Auffang- 
ebene Z= const. sia < X <0,—w < Y<o, 
erhalten wir : 


X, Y, Z sind dabei die Koordinaten des Bildraumes, 
J(R, Z) ist die Intensititsverteilung im Beugungsbild 
eines Objektpunktes in der Auffangebene Z = const., 
R der Abstand vom Zentrum, X, ist die X-Koordi- 
nate des Aufpunktes, d. h. sein Abstand von der 
Kante. Alle Gréssen sind wie in [2] normiert, auch 
1/4x ist ein Normierungsfaktor. An eine Integration 
von (1) auf analytischem Wege, wie sie im aberra- 
tionsfreien Fall von W. WEINSTEIN [3] erreicht wor- 
den ist, kann bei dem hier in Betracht zu ziehenden 
« grossen » Offnungsfehler nicht gedacht werden. Es 
diirften nur numerische Verfahren in Frage kommen. 
E. WANDERSLEB [4] nimmt eine feine Unterteilung 
des Integrationsgebietes vor und berechnet (1) durch 
reine Summation zunichst iiber Y (Lichtlinie) und 
dann tiber X. Da er J nur geometrisch-optisch in 
Ansatz bringt, J also ausserhalb des geometrisch- 
optischen Zerstreuungskreises verschwindet, wird in 
(1) das Integrationsgebiet endlich. Es treten dann 
keine prinzipiellen Schwierigkeiten auf. Mit den glei- 
chen Voraussetzungen arbeitet auch T. pi Francra|5]. 
Die wellenoptische Struktur der Lichtverteilung 
beriicksichtigt E. Fryxe [6]. Er entnimmt J experi- 
mentellen Messungen, approximiert stiickweise durch 
Polynome unter Beschrinkung auf die inneren vier 
Ringe, der dussere, experimentell auch schwer zu 
erfassende Lichtschleier wird vernachlissigt, und inte- 
griert dann (1) auf analytischem Wege. Wie [ 2 ], Abbil- 
dung 7 zeigt, enthilt aber bei grossem Offnungsfehler 
der dussere Lichtschleier noch einen betrachtlichen 
Energieanteil, so dass wir ihn hier voll beriicksichti- 
gen wollen. 


18 a 


Transformieren wir zuniichst (1) auf Polarkoordina- 
ten 


{ Aa | 5 Po 
7 


Ql) 2k = | ay J(R,Z) RAR de = 
2 ; ahaa pew istry. ela 
Pre (| Xeleos | re 
| | 7 J(R,Z)RdR ve]. te 1Ur 2e oS 0 
e 9 e/ 0 


und fiihren damit die Berechnung von «J, auf die 
Totalbeleuchtung 


id 


(3) BIR, Z)=— | SR", Z) RB! aR’ 
0 


a 
e 


zuriick, was fiir unsere Methode charakteristisch ist. 
B gibt die Totalbeleuchtung in einem Kreis mit 
Radius R in der Auffangebene Z = const. an und 
ist entsprechend der Energiebilanz normiert auf 

(4) B, = Bin, Z) = H2 =1, 


R 


H, = 1: normierte Einfallshohe des Randstrahls. 


Damit geht (2) iiber in 


{ T/2 We 
(5) J,(X,,Z) = — E +B 


mT Jo 


Nun zur Berechnung von B(R,Z). In [2] wurden 
fiir die Intensitatsverteilung J sowohl fiir achsennahe 
wie achsenferne Aufpunkte Formelsitze hergeleitet 
und eine Anzahl Beispiele numerisch durchgerechnet. 
Wir wollen zunachst als Ergiinzung hierzu nach 
N. G. van Kampen [7] noch einen Niherungsaus- 
druck aufstellen, welcher die Intensitat in der lichtsei- 
tigen Umgebung der Schattengrenze (Randstrahl) gut 
wiedergibt (Randstrahl selbst ausserhalb der Brenn- 
flachen). 


~ if 4 \2 
GRAIN pent | CH: 
| Pax Prr| \ 
(1 e 
A oh a= (8) et 
2 rei 
Alle Bezeichnungen wie in [2], Abschnitt 5. C(t) 


und S(t) sind die FResnetschen Integrale. Damit ist 
J(R, Z) in den wichtigsten Querschnitten Z = const. 
bis zur geometrischen Schattengrenze vollstindig 
bekannt und wir konnen B(R,Z) entsprechend (3) 
leicht durch numerische Integration (Simpsonsche 
Regel) berechnen. Im geometrischen Schatten kénnen 
wir durch riickliufige Integration von R = aus fiir 
B, — B(R,Z) genau in der gleichen Weise, wie es 
vom Verfasser schon im aberrationsfreien Falle durch- 


FOCKE 
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gefiihrt worden ist [8], direkte asymptotische Niihe- 


rungsausdriicke aufstellen. 
set a 
wv Pr 
— Hi), t=+4/——, 
Tv " 


Orr 


(7) 1—B(R,Z) ~ 


fiir R in der Ausseren Umgebung der Schattengrenze 
(diese selbst ausserhalb der Brennfliche), Hilfsfunk- 
tion H(t) entsprechend [8], Formel (11) und Tabelle 1. 


(8) t{—BR2) oo 


fiir R im Schattengebiet. R = R ist der Abstand des 
Randstrahldurchstosspunktes von der <Achse. Im 
Falle sphirischer Randkorrektion ist R = 2 Z, sonst 
folgt R unmittelbar aus [2], Formel (5, 2). 

Wir kommen somit bei der Berechnung von B von 
zwei Seiten an die Schattengrenze heran, und der 
Vergleich der beiden erhaltenen Werte gibt eine gute 
Kontrolle aller vorangegangenen Rechnungen. In 
den von uns durchgerechneten Beispielen bheb die 
Differenz unter 0,018, . Im Gegensatz zu der stark 
oszilherenden Intensititsverteilung haben wir fiir die 
Totalbeleuchtung eine wesentlich « glattere » Werte- 
verteilung erhalten. Diese gestattet nun leicht die 
Integration (5) wiederum numerisch auszufiihren. 
Die Stmpsonsche Regel liefert schon bei einer Unter- 


5 fire \e amas 
teilung des Integrationsintervalles (9; > )in3 x 2Teil- 


intervalle ausreichende Genauigkeit. Fiir grosse xX, 
kann man durch Einsetzen von (8) das Integral (5) 
auch geschlossen ausfiihren und erhalt asymptotisch 


(9) F(X 2) ~ = 


| fae 
xX arctg \/ —— 
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Im Falle der aberrationsfreien Abbildung geht (9) in 
die schon von WEINSTEIN angegebene Entwicklung [ 3] 
Forme! (19) iiber. 

Ist die Lichtverteilung J x(X,, Z) im Bild einer Halb- 
ebene bekannt, so kann man durch weitere Superpo- 
sition die Abbildung anderer inkohiarenter Objekt- 
strukturen untersuchen. Wir betrachten zunichst einen 
einzelnen hellen Streifen auf schwarzen Hintergrund. 
Sein geometrisches Bild sei— 3< NX <3, —w < Yeu 
mit der Intensitét 1, dann erhalten wir fiir die Inten- 
sitatsverteilung in seinem Beugungsbild 


(10) Is.(XJ=t—JI,(X, +8) —T,(8—X,) = 
= J,(X.—8) —I,(X, + 8), 
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Ags. 1. — Die Lichtverteilung im Bild einer Kante und einzelner Streifen verschiedener Breite unter 


dem Einfluss des Offnungsfehlers. (1): Das geometrische Bild der Testobjekte in der Bildebene. (2)-(5) : sukzessive 
anwachsender Offnungsfelher bis zur Wellenaberration Pimnax = 8 A, randstrahlkorrigiert. (6)-(8) : Offnungstehler 
Vinax = 84, jeweils um einen Blendenwert abgeblendet, pg : relativer Blendenradius, Objektintensitat mit e 
ansteigend, so dass in allen Fallen das gleiche eaitridehe Bild (1). Die Auffangebene liegt in allen Fallen 
im jeweiligen Hauptmaximum der Punktkaustik, [2], Abb. 1 und Abb. 6, Abszissen in normierten Bilde- 
beneneinheiten, 1 entspricht 2//27h, mit 2h,/f : Offnungsverh4ltnis (bei voller Offnung), A: Wellenlange 
des Lichts. 
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X , ist wieder die X-Koordinate des Aufpunktes von 
der Mitte des Streifens aus gerechnet, die Abhingig- 
keit von Z ist im Argument unterdriickt. Durch unend- 
lichfache Uberlagerung kénnen wir 
Abbildung eines Streifengitters iibergehen, 
breite wie oben 2 8, Streifenabstand 2 8, 
lange 1 = 48, 


weiterhin zur 
Streifen- 
Perioden- 


(11) Igy (Xy + 49 8). 


Die Aufpunktkoordinate ziihlt wieder von der Mitte 
des nullten Streifen aus. Fiir den wie iiblich durch 


J tes 


max min 


erklirten Gitterkontrast folgt dann 


c 


“max min 


(12) C,=C,(N)=1—4.J,(8) +4 J, (38)—4.J,(5 3)+... 

abhangig von 3 bzw. von der Linienfreqenz N = 1/48 

des Gitters. Fiir grosse § kann man entsprechend (9) 
2 1 

niherungsweise J,.(v5) ~ — — einsetzen und sum- 
qT v 

mieren. 


™O 
Schliesslich kénnen wir hieraus nach J. \W. Co.r- 
MAN [9] auch zur eigentlichen Kontrastiibertragung 
tibergehen, worunter wir nach P. LinpBere [1] den 
Bildkontrast eines Gitters mit sinusférmiger Hellig- 
keitsverteilung und Kontrast 1 verstehen wollen. 


= 3 [ow 


o-- 


rl 
(14) + 3 €(8N)— = C.(5N) + 


C(4.N) —...| 


Wie schon oben gesagt charakterisiert die Kontrast- 
tbertragung die Abbildung eines solchen Gitters 
vollstandig. Ausserdem gestattet sie durch Fourierent- 
wicklung zu anderen periodischen und durch Fourier- 
transformation zu nichtperiodischen linearen Objekt- 
strukturen tiberzugehen. 


C,(7 N) — Tre 
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Zum Abschluss der Arbeit wollen wir einige prak- 
tische Resultate mitteilen. Abbildung 1 zeigt die 
wellenoptisch berechnete Lichtverteilung im Bild 
einer Kante und einzelner Streifen verschiedener 
Breite unter dem Einfluss des Offnungsfehlers. Dieser 
bewirkt mit anwachsender spharischer Zone bzw. 
Wellenaberration einmal eine Abflachung des Licht- 
abfalls an der Kante selbst zum anderen eine bedeu- 
tende Zunahme des Licht- bzw. Schattensaumes zu 
beiden Seiten. Dabei ist der letztere Effekt stark 
dominierend, es geniigt deshalb keinesfalls, die Licht- 
verteilung allein durch den Intensititsgradienten an 
der Kante selbst zu charakterisieren. Dieser liisst sich 
unmittelbar aus (1) durch Substitution X’ = Y — NX 


Q 
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[Orr AcTa 
und Dilferentiation nach der neuen oberen Grenze 
— X, berechnen. 

Ce | i 
(15) dX, ef (O) ye | J(R) dk. 


Die entstehende Kontrastminderung wird durch die 
Betrachtung der Streifenabbildung noch deutlicher. 
Wahrend bei den breiten Streifen mit wachsendem 
Offnungsfehler zuniichst eine Verwaschung der Strei- 
fenkanten eintritt, nimmt nach kleinen Streifenbrei- 
ten zu auch die Kontrastwiedergabe betrichtlich ab. 
Die letzten drei Abbildungszeilen zeigen die Wirkung 


der Abblendung. Diese driickt zunichst den Licht- 
bzw. Schattensaum wesentlich herab, auch der Licht- 


abfall an der Kante selbst wird etwas steiler, bis der 
Einfluss des Offnungsfehlers ganz verschwunden ist, 
sich dann aber der Beugungseffekt bildverschlech- 
ternd bemerkbar macht. Die giinstigste Blendenwir- 
kung ist wohl zwischen pe? = 1/4 und 1/8 zu suchen, 


Abbildung 2 stellt die Kantenabbildung in verschie- 
denen Auffangebenen dar. Den steilsten Lichtabfall 


0 oe T r a — 1 1 
100 50 ; 
Ass. 2, — Die Lichtverteilung im Bild einer Kante bei 


randkorrigiertem Offnungsfehler ©)... = 8 Ain den Auffang- 


ebenen Z=0,75A bisiA.Z=A : Zonenstrahlschnittpunkt, 
: Grésse der pee Zone, Z-Achse in Lichtrichtung, 
Abszissen wie Abb, 


finden wir in Ubereinstimmung mit dem Hauptmaxi- 
mum der Punktabbildung be Z = 0,9 A. Dagegen 
zeigt sich bei Z = 0,8 A ein Riickgang des diusseren 
Lichtsaumes. Da beide Effekte die Bildgiite beein- 
flussen, muss die giinstigste Einstellebene wohl zwi- 
schen beiden Z-Werten liegen. 

Schliesslich zeigt noch Abbildung 3 den Einfluss 
des Offnungsfehlers an Hand der Kontrastiibertra- 
gung. F ordert man z. B. eine Kontrastwiedergabe 
C = 0,5, so sinkt das Auflésungsvermégen mit 
Sachin Offnungstehler von 1! = _ im aberra- 
tionsfreien Falle auf 1 =70 bei © 8 2. Die Kon- 


my ax 
trastiibertragung verbindet die beiden fiir die Bildgiite 
ausschlaggebenden Komponenten Kontrast und Auflé- 
sungsvermogen zu einer sehr instruktiven zahlenmiissi- 
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up gen Leistungsangabe. Nach den kleinen Gitterperio- 
den zu treten allerdings Fehlangaben, sog. falscher 
Kontrast auf (in Form von Kontrastminima und 
-maxima), welcher nur durch die Gitterstruktur des 
085A Testobjektes hervorgerufen wird und keine allge- 
084 meine Aussage tiber die Bildgiite darstellt. Dieser 
O74 Effekt zeigt sich auch an den von P. LinDBERG expe- 
rimentell gemessenen Kontrastkurven [1] Figure 12, 
welche im iibrigen mir unseren theoretisch berechneten 
prinzipiell gut tibereinstimmen. 


LITERATURHINWEISE 


[1] P. LinpBere, Opl. Acta, 1, 1954, p. 80-89. 
[2] J. Focxe, Opt. Acta, 8, 1956, p. 110-126. 
A HAY (3] W. WEINSTEIN, J. O. S. A., 44, 1954, p. 610-615. 
[4] E. WANpDERSLEB, Die Lichtverteilung in der axialen Kau- 


x stik eines mit sphdrischer Aberration behafteten 
Objektivs, Berlin 1952. 


z f ° y al . 5] T. pi Franera, J. O. S. A., 48, 1953, p. 827-835. 
ABB. 3. — Die Kontrastitbertragung C in Abhiangigkeit tet E. Finke, Jenaer Jahrbuch 1954, IT. Toil, p. 385-400. 
vom Offnungsfehler. Entspricht den Fallen (2)-(5) von Abb. 1. [7] N. G. vAN Kampen, Physica, 14, 1948, p. 575-589. 
l: Periodenlange des Gitters in normierten Bildebeneneinhei- [8] J. Focke, Opt. Acta, 3, 1956, p. 161-163. 
ten wie in Abb. 1. [9] J. W. Cottmann, J. O. S. A., 44, 1954, p. 468-471. 


Ein Verfahren zur Messung der Kontrast- 
Ubertragungsfunktion photographischer Objektive. 


K. Rosenuaver und K.-J. Rosensprucu, Braunschweig 
Mitteilung aus der Physikalisch-Technischen Bundesanstalt 


Summary. — The contrast Kx of a self-luminous or incoherently illuminated object is defined by K = (Jmax — Jmin)/(Jmax+ Jmin) 
where J is the intensity of illumination. When an optical system forms an image of the object, the contrast is reduced, due to 
diffraction, aberrations, aid veiling-glare. The reduction of contrast is a function of the size of the object (or the spacial frequency), 
the relative aperture of the imaging system and its state 6f correction. 

An apparatus is described for measuring this contrast transmission function for miniature camera lenses. The objective to be 
tested projects on to a narrow slit the image of a rotating grating having varying frequencies, and K = 1. A photomultiplier 
gives an output proportional to the intensity of transmitted light and a cathode ray oscillograph synchronised with the rotation 
of the grating gives a stationary representation of the light distribution in the tmage. 

Results obtained with this apparatus from miniature camera lenses are given and discussed for radial and tangential lines at 
various relative apertures and field angles and at different axial planes of focus. 


SomMatreE. — Le contraste K d’un objet éclairé de facon incohérente ou lumineux par lui-méme est défini, si, J désigne Vintensité 
par K = (Jmax — Jmin)/(Jmax + Jmin). Lorsqu’'un systéme optique forme une tmage de cet objet, le contraste subit une dimi- 
nution due a la diffraction par Vouverture, aux aberrations de optique et a la lumiére parasite diffusée. Cette baisse de contraste 
est une fonction de la grandeur de lV objet (ou de la fréquence spatiale) et dépend de Vouverture relative et de V’état de correction du 
systéme optique. 

Un dispositif est élaboré pour la mesure de cette fonction de transmission du contraste dans le cas des objectifs photographiques 
de petit format. Avec Vobjectif a essayer, on projette sur une fente étroite Vimage d'un réseau tournant, de pas décroissant et de 
contraste K =1. Un photomultiplicateur convertit la lumiére transmise en tensions proportionnelles aux intensités, et sur un oscillo- 
graphe cathodique synchronisé avec la rotation du test on oblient une image fixe dela répartition des intensités dans V image du test. 

On présente et on discute les courbes de transmission de contraste obtenues avec ce dispositif pour des objectifs photo de petit 
format, dans le cas de réseaux 4 traits horizontaux et verticaux pour diverses valeurs de Pangle et du diaphragme, ainsi que dans 
le cas d’un défaut de mise au point dans Vaze. 

ZUSAMMENFASSUNG. — Der Kontrast K eines inkohdrent beleuchteten oder selbstleuchtenden Objektes, definiert nach der Gleichung 
K = (Jmax — Jmin) | (Imax + Jmin) (J = Intensitdi), wird bei der Abbildung durch optische Systeme durch Beugung an der 
Offnung, durch die Abbildungsfehler der Optik und durch nicht bilderzeugendes Streulicht verringert. Diese Abnahme des 
Kontrastes ist eine Funktion von der Grdésse der abgebildeten Objekte (Raumfrequenz) und hdngt von dem Offrungs 
verhdltnis und dem Korrektionszustand des abbildenden Sysiems ab. 

Fir die Messung dieser Kontras!-Ubertragungsfunktion wurde fiir Kleinbildobjektive eine Apparatur entwickelt. Ein 
rotierendes Gitler mit abnehmender Gitterkonsltante und einem Kontrast K = 1 wird durch das zu priifende Objektiv auf 
einen engen Spalt abgebildet. Das durchgelassene Licht wird von einem Photo-Muliplier in intensitdtsproportionale Span- 
nungen verwandelt, die auf einem mit dem rotierenden Test synchronisierlen Kathodenstrahl-Oszillographen eine stehendes 
Bild der Intensitdtsverteilung des Testbildes erzeugen, Es werden die mit dieser Apparatur fiir Kleinbildobjektive erhaltenen 
Kontrastitbertragungskurven fiir horizontale und vertikale Gitterlinien bei verschiedenen Bildwinkeln und Abblendungen, sowie 
fiir Defokussierung in der Achse gezeigt und diskutiert, 


Bei der Priifung photographischer Objektive unter- denen das Zusammenwirken der Einzelfehler auf die 
scheidet man grundsitzlich zwei Methoden. Gesamtqualitiit des Bildes ermittelt wird. Das Ver- 
Erstens die differenzierenden Methoden, bei denen fahren, das hierbei hauptsichlich angewandt wird, 
die einzelnen Bildfehler der Objektive bestimmt wer- ist das Verfahren der Beurteilung der Bildqualitat 
den und zweitens die integrierenden Methoden, bei einer Testtafel durch visuelle Beobachtung des Luft- 


22 K. 


oder Mattscheibenbildes oder durch Beurteilung einer 
photographischen Probeaufnahme. Die Testtafel ist 
aus geeigneten Testen verschiedener Grosse und ver- 
schiedenen Kontrastes aufgebaut. Als Kriterium fiir 
eine gute Bildqualitat sieht man meistens das Auflo- 
sungsvermégen an. Unter Auflésungsvermégen wird 
dabei der Abstand zweier Punkte verstanden, die 
gerade noch getrennt wahrnehmbar sind. Bei einem 
gitterartigen Test wird analog dazu die Gitterkon- 
stante oder auch ihr reziproker Wert als Auflosungs- 
vermégen bezeichnet. Es ist schon lange aufgefallen, 
dass bei manchen Objektiven, obgleich ein hohes 
Auflésungsvermégen erzielt wurde, das Gesamtbild 


doch merkwiirdig kontrastarm oder — wie gesagt 
wird — das Gesamtbild wattig oder weich erscheint. 


Es drangt sich so die Frage nach dem Zusammen- 
hang zwischen Auflosungsvermogen und Kontrast 
auf. Setzt man voraus, dass als Objekt z. B. ein Gitter 
aus vollkommen undurchsichtigen Strichen und voll- 
kommen durchsichtigen Zwischenriumen verwandt 
wird, so ist der Kontrast gleich 100 °%, wenn man fol- 
gende Kontrastdefinition verwendet : 


Max” J Min 


100 x 
I oe Min 


aa KO OY. 
K%. 


Man wird dann bei abnehmender Gitterkonstante 
Auflésung bis zu dem Punkt erzielen, wo man die 
Kontrastschwelle des Auges fiir den vorliegenden 
Beobachtungsfall erreicht. Die Kontrastschwelle des 
Auges ist aber von vielen Nebenumstanden, wie Ver- 
grosserung, Beleuchtungsstarke, Giite der Trennkante, 
Farbe und Blendung abhangig. Man wird also nicht 
erwarten kénnen, dass Auflésungsvermégen photo- 
graphischer Objektive, die von verschiedenen Beob- 
achtern nach verschiedenen Verfahren  ermittelt 
wurden, iibereinstimmen. Es sei denn, man gibt die 
bei der Messung wirksamen Randbedingungen genau- 
estens an. 

_ In der elektrischen Nachrichtentechnik ist die 
Ubertragungstheorie eines Vierpols weitgehend aus- 
gebildet worden und man kann fiir den Fall der opti- 
schen Ubertragung eines Gitters mit abnehmender 
Gitterkonstanten ahnliche Betrachtungen wie im elek- 
trischen Fall anschliessen [1]. Auch in der Optik 
werden von einem Sender (Objekt) Frequenzen (in 
der Optik verschiedene Gitterkonstanten auch als 
Raumfrequenzen bezeichnet) mit verschiedenen Amphi- 
tuden und Phasen ausgestrahlt, die durch ein Linsen- 
system (Ubertrager) nach Phase und Amplitude ver- 
indert werden und am Empfangsort (Bildebene) 
durch Uberlagerung das Abbild des Senders erzeugen. 
Ks interessiert nun, mit welchen Amplituden und 
manchmal auch Phasen die verschiedenen Raumfre- 
quenzen iibertragen werden. Ein Mass dafiir ist das 
Verhaltnis des Kontrastes im Bild zu dem im Objekt 
bei verschiedenen Raumfrequenzen. In der Nachrich- 
tentechnik wiirde man _ vielleicht vom Modulations- 
grad sprechen. Hier wird dieser Zusammenhang als 
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Ubertragungsfunktion bezeichnet. Diese Ubertra- 
gungstheorie ist unter anderem von Linroot [2], 
H. H. Hopkins [3] und Mitarbeitern [ 4] ausgearbei- 
tet worden. 

Grundsiitzlich kann man fiir die Ubermittlung der 
Ubertragungsfunktion zwei Wege beschreiten. Erstens 
kann man auf der Eingangsseite Frequenzen konstan- 
ter Amplitude eingeben und kann auf der Ausgangs- 
seite den Zusammenhang zwischen Amplitude und 
Frequenz ermitteln. Zweitens ist es moglich, die Ein- 
gangsamplitude bei verschiedenen Frequenzen so zu 
regeln, dass die Ausgangsamplitude einen konstanten 
Wert annimmt. Mit beiden Methoden kann man die 
Ubertragungsfunktion nach Amplitude und Phase in 
Abhingigkeit von der Frequenz ermitteln. Diese Ver- 
haltnisse sind der Nachrichtentechnik weitgehend 
bekannt. Das in der Optik iibliche Verfahren zur Be- 
stimmung des Auflésungsvermégens arbeitet mit einer 
konstanten Ausgangsamplitude, naimlich der Kontrast- 
empfindlichkeit des Auges in den besonderen Beobach- 
tungsfillen. Man kann also fiir Variationen des EKingang- 
kontrastes bei Beobachtung verschieden weiter Gitter 
so die Ubertragungsfunktion, oder in der Optik 
richtiger Kontrast-Ubertragungsfunktion genannt, 
ermitteln. Da aber das Auge als _ physiologisches 
Messinstrument starken Schwankungen sowohl von 
Beobachter zu Beobachter wie auch zeitlichen Schwan- 
kungen usw. unterliegt, so ist es klar, dass solche 
Messungen mit relativ grossen Fehlern behaftet sind. 


Man hat deshalb nach verbesserten Verfahren 
gesucht, und es wird in diesem Zusammenhang auf 
die Arbeiten von INGELSTAM u. LINDBERG [5], 
ScHADE [6] [7] [8], SELvyn u. TEARLE [9], Linp- 
BERG [10] und H. H. Hopkins [11 ] hingewiesen. An 
ein fiir die Praxis verwendbares Verfahren wird 
aber die Bedingung gekniipft, dass es neben ausrei- 
chender Genauigkeit auch geniigend schnell durch- 
fiihrbar ist und dass die Apparatur mechanisch nicht 
zu schwierig herzustellen bzw. nicht zu kostspielig ist. 
Ausserdem miissen die Ergebnisse schnell greifbar 
sein. 

Es wurde deshalb folgender Weg beschritten, der in 
Abb. 4 schematisch dargestellt ist. In der Objektebene 
eines Photoobjektivs, die ca 5 m von dem zu priifen- 
den optischen System entfernt war, befand sich ein 
Gitter, dessen Linien scharfe Kanten zeigten und des- 
sen Kontrast gleich 100° bzw. 1 war. Wiirde man 
nun in der Bildebene die Lichtverteilung z. B. mit 
Gleichlicht ermitteln, so kénnte man das Verhiltnis 
zwischen dem Kontrast des Objektgitters und dem 
Kontrast seines Bildes ermitteln. Da dieses Gitter eine 
in geometrischer Reihe abnehmende Gitterkonstante 
besitzt, so wiirde man daraus also den Zusammenhang 
zwischen Kontrastibertragung und Auflésungsvermé- 
gen ermitteln kénnen. Eine solche Anordnung hat aber 
den Nachteil, dass waihrend der gesamten Messzeit 
sich die Helligkeit der Beleuchtungslampe nicht 
andern und dass auch die Empfindlichkeit des zur 
Registrierung benutzten Empfangers keinen Schwan- 
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rotierendes 
Gittertest 


kungen unterliegen darf. Als weitere Schwierigkeit 
kommt hinzu, dass man die Abtastung des in der Bild- 
ebene entstandenen Luftbildes mit einem sehr engen 
Spalt vornehmen inuss. Das Auflésungsvermégen von 
handelsiiblichen photographischen Kleinbildobjekti- 
ven von ca 50mm Brennweite kann aber bis zu Werten 
von etwa 4 » gehen. Man miisste also als Abtastspalt 
einen solchen von wenigstens 0,8 bis 1 uw verwenden 
und diesen Spalt mit der entsprechenden Genauigkeit 
an dem in der Bildebene befindlichen Gitterbild ent- 
lang fiihren. Es wiirden also sehr hohe Anforderungen 
an die Giite einer solchen Fithrung und an die Stabili- 
tat der gesamten Anordnung gestellt werden. Benutzt 
wurde deshalb ein Spalt von tblichen Grossenabmes- 
sungen, in diesem Fall von etwa 3/100 mm, der durch 
ein gutes, ausgesuchtes Mikroobjektiv etwa 40 fach 
verkleinert in der Bildebene abgebildet wird. Dieses 
Mikroobjektiv (num. Apertur 0,7) bildet infolge seines 
grossen Offnungswinkels auch die unter grossen Bild- 
winkeln einfallenden Biindel auf dem Spalt ohne Dre- 
hung des Mikroskops ab, wenn das Test auf einer die 
Gegenstandsebene darstellenden optischen Bank seit- 
lich verschoben wird. Natiirlich ist danneine seitliche 
Verschiebung des Mikroskops mit Spalt und Multi- 
plier notig. Abb. 2 gibt eine Ansicht dieser Einrich- 
tung, auf der das zu priifende Photoobjektiv, das 
Mikroobjektiv, der Photomultiplier, die seitliche wie 
auch die axialen Verschiebungsmoglichkeiten zu 
sehen sind. Die vorn sichtbare Messuhr dient zur Be- 
stimmung der axialen Verschiebung. Anstelle einer 
Bewegung des Spaltes wird eine Bewegung in der 
Objektebene vorgezogen. Es wird also das Gitter in der 
Objektebene bewegt und man erhilt auf dem hinter 
dem Spalt befindlichen Photomultiplier einen Strom, 
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Oszillograph 


Ass. 1. — Schema der Versuchsanordnung. 


2. — Prifeinrichtung. 


ABB. 


der der Lichtverteilung im Bild entspricht. Das ver- 
wandte Mikroobjektiv ist in seiner Qualitiét — wie 
sich gezeigt hat — wesentlich besser als die zu unter- 
suchenden photographischen Objektive. Im allgemei- 
nen kann man sagen, dass die geometrisch optischen 
Bildfehler guter Mikroobjektive etwa 10 mal so klein 
als die der hier zu untersuchenden photographischen 
Objektive sind. Sein Einfluss wird daher vernach- 
lassigt. Genauere Untersuchungen laufen noch. Um 
nun von den Schwankungen der Lichtquelle und von 
der Instabilitét moglichst unabhingig zu sein, wurde 
das Test auf einer schnell umlaufenden Trommel (ca 
20 sec!) (Abb. 3) angebracht. Man erhalt dann auf 
einem hinter dem Photomultiplier befindlichen Oszil- 


lographen unter Zwischenschalten geeigneter Ver- 
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ABB. 3. — Testtrommel mit Beleuchtung auf Schiene 
verschiebbar, drehbar und um 90° schwenkbar. 


starker ein stehendes Bild, wenn man die Zeitablenkung 
des Oszillographen mit dem Trommelumlauf synchro- 
nisiert. Man sieht also auf dem Oszillographenschirm 
die Lichtverteilung in der Bildebene des photogra- 
phischen Objektivs (Abb. 4a und 4b). Die rotierende 
Trommel ist um die Bildmitte um 90° drehbar, so dass 
sowohl horizontale als auch vertikale Gitterlinien 
abgebildet werden kénnen. Verstairker und Oszillo- 
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graph miissen so gewaihlt werden, dass sie die fur 
die Ubertragung notwendigen Frequenzen einwand- 
frei wiedergeben, was mit Rechteckimpulsen geeigne- 
ter Frequenz gepriift worden war. Die Beleuchtung 
des Testes erfolgte im durchgehenden Licht mit aus- 
reichender Apertur und zwischengeschalteter Streu- 
scheibe, um die Inkoharenz des Lichtes zu gewahrleisten 
und um nicht scharfe Interferenzstrahlen, die bei Git- 
tertesten eine fehlerhafte Messung des Auflésungsver- 
mdégens hervorrufen, zu erhalten. 

Die vorliegenden theoretischen Untersuchungen [ 3 } 
[4 | beziehen sich fast alle auf Teste mit einer Intensi- 
tatsverteilung, die sinusférmig und unendlich lang ist. 
Praktisch lassen sich solche sinusfoérmigen Verteilun- 
gen aber nur schwierig herstellen, man zieht deshalb 
Gitter vor, die rechteckige Lichtimpulse erzeugen. In 
diesem Fall muss man ein Gitter mit verschiedenen 
Gitterkonstanten benutzen. Hier wurden rechteckige 
Lichtimpulse benutzt, die durch einen Test erzeugt 
wurden, dessen Gitterkonstanten in einer geometrischen 
Reihe abnahmen. Jeweils ein schwarzer und ein 
weisser Strich hatten dabei gleiche Breite. Im Prinzip 
kénnte man natiirlich auch als Testobjekt Gruppen 
gleicher Gitterkonstante verwenden, die hintereinander 
angeordnet sind. Grundsatzlich ist auch die Verwen- 
dung einer einzelnen Kante méglich, denn eine solche 
scharfe Kante stellt ein Fourierspektrum dar, in dem 
alle Frequenzen vorkommen. Allerdings sind die héheren 
Frequenzen dabei nur mit kleinen Amplituden vertreten 
und deshalb ist eine solche Kante als Objekt schlecht 
geeignet. Da das Fourierspektrum fiir die hier verwen- 
deten Teste zur Zeit noch untersucht wird und der Ober- 
wellengehalt der einzelnen Gitterstriche nicht genau 
bekannt ist, so wird auf eine Bewertung der Kanten- 
steilheit bei der hier benutzten Methode nur geringer 
Wert gelegt. Es muss noch bemerkt werden, dass 
beim Vergleich mit den von H. H. Hopkins [3] 
errechneten theoretischen Kurven zu beachten ist, 
dass bei den hier durchgefiihrten Messungen weisses 
Licht benutzt wurde. Die als Lichtquelle benutzte 
Gliihlampe hatte eine Farbtemperatur, die zwischen 
2 000 °K und ca 3 000 °K gehalten wurde. In Verbin- 


Ass. 4. — Bild auf der Kathodenstrahlréhre der yom Objektiv tibertragenen 
Intensitatsamplituden des Testes. z 


a) gut korrigiertes Objektiv Blende 2,8, 
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b) einfacheres Objektiv Blende 2,8. 


Linke Seite : Bild des gesamten Testes ; Rechte Seite: Bild des Testes yom 10. 
Gitterstrich ab. (Zeitachse ist elektrisch gedehnt). 


dung mit der spektralen Empfindlichkeitskurve des 
benutzten Photomultipliers ergab sich so eine spek- 
trale Empfindlichkeit, die etwa der entsprach, die man 
fiir panchromatisches Photomaterial in Kombination 
mit einem schwachen Gelbfilter fiir normales Tages- 
licht erhalt. Beim Vergleich der theoretischen Kurven 
mit den hier erhaltenen, kann man immer nur eine 
annihernde Ubereinstimmung erwarten, denn, da die 
Objektive eine chromatische Liingsaberration besitzen, 
so iiberlagern sich verschiedene Ubertragungskurven 
und es kann immer nur fiir eine Farbe eine wirkliche 
Fokussierung erfolgen. Fiir die davon abweichenden 
Farben erhalt man dann die Ubertragungsfunktion, 
wie sie sich fiir eine defokussierte Einstellung ergibt. 

Die Ubertragungsfunktion selbst wurde so gebildet, 
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dass die auf dem Oszillographenschirm  sichtbare 
Kurve photographiert wurde und in der Projektion 
nach der vorne angegebenen Formel der Kontrast als 
Funktion des Gitterabstandes ermittelt wurde. Um 
die Messungen der Praxis anzugleichen, und um einen 
einheitlichen Bezugspunkt zu haben, wurde dabei die 
Einstellung auf bestes Auflésungsvermégen fiir das 
grésste Offnungsverhaltnis auch fiir die anderen 
Offnungsverhiiltnisse, wie auch fiir die von der Achse 
abweichenden Bildwinkel beibehalten. Man kann so 
durch Ablesung an der in Abb. 2 sichtbaren Messuhr 
durch Nachstellen die Blendendifferenz ermitteln und 
auch messen, wie die Bildfeldwélbung verlaiuft. Daman 
diese Messungen sowohl fiir den meridionalen wie auch 
fiir den sagittalen Schnitt durch Schwenken der Test- 
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Apps. 5. — Kontrastitbertragungsfunktion fiir ein IKleinbildobjektiv f= 50 mm, 1:2 ; 
Blende 2 senkrechte Gitterlinien fiir Bildwinkel von 0°, 8°, 16°, 24°. 


26 K. ROSENHAUER UND K.-J. ROSENBRUCH [Opr. Acva 
40— 
AS 
os: 
. aN ath 
N \ ane Bildwinkel 
o7. Hy \ : 


06 


—e 
S 
a 


02 


Kontrastabertragungsfaktor 
S 


of 
4 | | inte \ Sees | See | 
0 50 100 150 200 250 300 350 
Raumfrequenz in Linien/mm 
ABB. 6. — Kontrastiibertragungsfunktion fir ein Kleinbildobjektiv f = 50 mm, 1: 2; 


Blende 4 senkrechte Gitterlinien fir Bildwinkel von 0°, 8°, 16°, 24°. 


trommel und des Spaltes vor dem Multiplier um einen 
Winkel von 90° durchfiithren kann, so kann man hier- 
mit auch den Astigmatismus feststellen. Abb. 5 stellt 
die Ubertragungsfunktion eines hochqualifizierten 
Kleinbildobjektives von 50 mm Brennweite und fiir ein 
Offnungsverhiltnis von 1:2 fiir die Achse und fiir 
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verschiedene Bildwinkel dar. Abh. 6 zeigt das gleiche 
fiir Blende 4,0. Abb. 7 zeigt die gleichen Verhaltnisse 
fiir den hierzu senkrechten Schnitt. Ausserhalb der 
Achse kann man einen Unterschied der Ubertragungs- 
funktion fiir den meridionalen und fiir den sagittalen 
Bildschnitt, d.h. den Einfluss des Astigmatismus 
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Blende 4 wangerechte Gitterlinien fiir Bildwinkel von 0°, 8°, 16°, 24°, 
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Azss. 8. — Kontrasttitbertragungsfunktion fiir ein Kleinbild- 


objektiv f = 50 mm, 1: 2; Blende 11 senkrechte Gitter- 
linien fir Bildwinkel von 0, 8°, 16°, 24°. 


bei Vergleich mit Abb. 6, erkennen. Abb. 8 zeigt die 
Ubertragungsfunktion fiir das gleiche Objektiv, jedoch 
fiir eine Abblendung auf 1:11. In Abb. 9 wird die 


Kontrastubertragungsfaktor 


a7 


Wirkung einer Defokussierung dargestellt. Die Kurve 0 
zeigt die Ubertragungsfunktion in der Achse fiir das 
in Abb. 5 erwihnte Kleinbildobjektiv und die mit —1, 
— 2 baw. + 1, + 2 usw. bezeichneten Kurven stellen 
die Ubertragungsfunktion fiir eine Defokussierung 
intra- bzw. extrafokal in 1/100 mm dar. 

Aus der Kurve der Ubertragungsfunktion z. B. aus 
Abb. 6 kann man — wie z. B. Abb. 10 zeigt — das 
Auflésungsvermégen des Objektivs bei der betreffen- 
den Abblendung fiir die Achse und fiir verschiedene 
Bildwinkel bei verschiedenen Kontrasten ermitteln. 
Aus der in Abb. 9 dargestellten Wirkung der Defokus- 
sierung liasst sich fiir eine vorgegebene Unschiarfe und 
fiir einen vorgegebenen Kontrast die Tiefenschirfe 
des Objektivs unter Beriicksichtigung seiner Korrek- 
tion ermitteln. Man kann so die mit Hilfe des Zer- 
streuungskreises fiir ein ideales Objektiv errechnete 
Tiefenschirfe mit der wirklich erreichten vergleichen. 

Weiterhin kann man die tatsichlich erreichte und 
in der Kontrast-Ubertragungsfunktion der Gesamt- 
wirkung sichtbar gemachte Abbildungsgiite mit der 
theoretisch fiir ein ideales Objektiv errechneten ver- 
gleichen. Allerdings ist hierbei zu beachten, dass die 
hier gemessene Kontrast-Ubertragungsfunktion fiir 
weisses Licht ermittelt wurde, wihrend der berech- 
neten eine monochromatische Lichtquelle der Wellen- 
lingen 0,5 » zugrunde gelegt wurde. Abb. 11 zeigt 
diese Verhiltnisse fiir ein hochqualifiziertes Objektiv 
(Kleinbildobjektiv f = 50 mm 1: 2). Die Abszisse ist 
hier jedoch nicht in Raumfrequenzen geteilt, sondern 


App. 9. — Kontrastiibertragungsfunktion eines Kleinbildobjektives / = 50 mm, 
1: 2 fiir die Achse und volle Offnung bei verschiedener Defokussierung. 
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ABB. 10. — Auflosung eines Kleinbildobjektives f = 50 mm, 1: 2 
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bei Blende 4 in Abhangigkeit vom Bildwinkel bei verschiedenem I<ontrast. 


es ist eine normierte Raumfrequenz zugrunde gelegt 
worden. Die normierte Raumfrequenz s leitet sich aus 
der Raumfrequenz R, der Apertur des Objektivs 
n sin x und der Wellenlainge % nach folgender Formel 
abs: 

mn 


n sin « 


Man sieht aus der Abb. 11, dass mindestens von der 
Blende 5,6 ab praktisch die theoretische Kurve 
erreicht wird. Abb. 12 zeigt das gleiche fiir ein einfa- 
cheres Photoobjektiv. Hierbei zeigt sich, dass in der 
Achse selbst die Blende 11 noch nicht die theoretische 
Kurve erreicht. Man miisste diese Verhiltnisse auch 
fiir verschiedene Bildwinkel untersuchen und kénnte 
so ein Gesamturteil iiber die Gesamtbildqualitiat eines 
Photoobjektives erhalten. 

Bedenkt man nun noch, dass die unter der theore- 
tischen Kurve befindliche Fliche der Definitionshel- 
ligkeit nach Srrent [12] entspricht, so stellt das 
Verhiltnis der Flicheninhalt> unter der gemessenen 


Kurve zur Gesamtfliche unter der theoretischen 
Kurve den praktisch ermittelten prozentualen Anteil 
der Definitionshelligkeit dar. 

Es mag noch erwahnt werden, dass in der prakti- 
schen Photographie auch die photographische Schicht 
eine Kontrast-Ubertragungsfunktion _ besitzt. Legt 
man die Werte zugrunde, die FrirsEr [ 13] ermittelt 
hat, so kann man daraus die Kontrast-Ubertragungs- 
funktion errechnen, die durch die gesamten Wirkun- 
gen von Objektiv und Photomaterial entstehen. Zu 
beachten ist dabei, dass bei Hintereinanderschalten 
zweier Ubertrager entsprechend den Ergebnissen der 
Vierpoltheorie eine Multiplikation der Kontrast-Uber- 
tragungsfunktionen erforderlich wird. Man sieht hier- 
durch leicht ein, dass eine lineare Superposition des 
Auflésungsvermégens der photographischen Schicht 
und die des Auflésungsvermégens des photographi- 
schen Objektives im allgemeinen nicht zu richtigen 
Werten fiihren kann. Abb. 13 stellt die Kontrast- 


Ubertragungsfunktion fiir das System aus einem gut 
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Apps. 11. — Normierte Kontrastittbertragungsfunktion fiir die Achse eines idealen Objektives 
und fir ein Kleinbildobjektiv f= 50 mm, 1: 2 bei verschiedener Abblendung.* 


—— ideales Objektiy 
Blende: 
--- 29 


as § 6 einfaches Kleinbildobjektiv 


14 


Kontrastubertragungsfaktor 


if ! 1 | 
0 02 Ot 06 08 
normierte Raumfrequenz 


Asp. 12. — Normierte Kontrastitbertragungsfunktion fiir die Achse eines idealen Objektives 
und fiir ein Kleinbildobjektiv geringerer Gitte f = 50 mm, 1: 2,9 bei verschiedener Abblendung.* 


* Bei einer Uberpriifung der Messwerte ergab sich, dass die in s eingesetzten Aperturwerte n- sin aw um 20 % 
7U gross waren, so dass bei Beriicksichtigung dieser KXorrektur alle ICuryven unter der idealen liegen. 
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Kleinbildobjektiv / = 50mm, 1: 2 bei Blende 4 und feinzeichnender Film. 


auflésenden handelsiiblichen Film und einem Wlein- 
hildobjektivy fiir die optische Achse dar. Bedenkt man, 
dass die Kontrastschwelle des Auges bei etwa 0,2 hegt, 
so sieht man, dass fiir das Gesamtauflosungsvermigen 
yon Objektiv und Film sich Werte von etwa 60-70 
Linien pro mm ergeben, was mit direkten praktischen 
Messungen gut iibereinstimmt. 


Aus Vorstehendem ergibt sich die praktische 
Brauchbarkeit der genannten Methoden. Zur Zeit 


wird diese Methode auch auf die Bestimmung des 
Auflisungsvermégens bzw. der Kontrast-Ubertra- 
gungsfunktion von Fernrohren angewandt. Uber die 
Ergebnisse wird spiter berichtet werden. 

Der Einfluss des nicht bilderzeugenden Lichtes, das 
durch diffuse Streuung im Objektiv z. B. durch Viel- 
fachreflexion, Triibheiten im Glas, Streuung an den 
Fassungsrindern u. i. entsteht, wird mit der angege- 
benen Methode nur in dem Masse erfasst, wie der Test 
selbst als Streulichtquelle wirkt. Unter ungiinstigen 
Umstiinden kann jedoch der ganze Halbraum vor dem 
Objektiv als Streulichtquelle wirken und das Bild mit 
einem Streulichtuntergrund versehen. Die gemessenen 
frequenzabhangigen Kontrastiibertragungskurven wer- 
den daher in der Praxis oft um einen frequenz- 
unabhingigen Betrag verindert, der vom Streufaktor 
des Objektivs abhingt und in seiner Grésse mit der 
Lichtverteilung zwischen abgebildetem Objekt und 
Streulicht erzeugender Umgebung stark schwanken 
kann. 

Ferner ist auch die Intensitiitswiedergabe in der 
Bildebene in Abhingigkeit yom Bildwinkel oft yon 
Interesse, die von der Transparenz und Vignettierung 


abhiingen. Die Kontrast-Ubertragungskurven der 
Abb. 5-8 geben jedoch nur das_ frequenzabhingige 
Verhiltnis der Kontraste bei der gemessenen Linien- 
frequenz zur Linienfrequenz Null an, wahrend das 
Bild 4 direkt die Intensitaten erkennen lisst. 

Zum Schluss méchten wir nicht verfehlen, unserem 
Mitarbeiter, Herrn Herm, fiir seine wertvolle Hilfe 
bei der Herstellung der Apparaturen wie auch bei der 
Aufnahme der Kurven unseren Dank zu sagen. Ebenso 
Herrn Dr Hopkins fiir die Anregungen, die sich aus 
den Diskussionen dieses Verfahrens ergeben haben. 
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EFFECT ON THE HUMAN ELECTRORETINOGRAM aL 


The effect on the human electroretinogram of the distribution 
of flux in a light stimulus of finite duration * 


Lucia Roncut and J. D. MoreLanp 


Istituto Nazionale di Ottica, Arcetri, Firenze (Italia) 


Summary. — The dependence of the human ERG on the flux distribution in a white light stimulus is examined for two observers in 


the dark adapted state. The stimulus is broken down into three variables : 


lime of rise, time of exposure, and peak intensity. The 


response is analysed with regard to three characteristics : height, slope and latency-lime of the b-wave. 
On increasing independently exposure time or peak intensity, both the height and slope of the b-wave increase and lend to satu- 
ration in the limit. Latency-lime shows little or no change. A BUNSEN-Roscor relation between peak-flux and exposure time 


appeared to hold for one observer. 


When the time of rise is changed, keeping exposure time and peak intensity constant, no change in the ERG is recorded although 
the total energy contained in the stimulus may be changed twofold. The analysis of the ERG with a view to revealing a BUNSEN- 
Roscor relation is therefore an oversimplifying procedure: the rod mechanism cannot be considered as a simple integrating system. 

Short-duration stimuli reveal the x-wave response (characteristic of the cone mechanism) which appears in three modes of rou- 


ghly equal frequency. 


Analysts of records obtained for standard stimuli over a period of months shows a typically seasonal variation of appreciable 


magnitude in the height of the b-wave response. 


SomMAIRE. — L’influence de la distribution du flux sur UVélectrorétinogramme (ERG) humain pour un stimulus de lumiére blanche 


est étudié pour 2 observateurs adaptés & UVobscurité. 


Les 3 composantes du stimulus sont le temps d’établissement, le temps d’exposition el le maximum d’ intensité. 

La réponse est analysée en fonction de trois caractéristiques : hauteur, pente et temps de latence de l’onde-b. 

Si le temps d’exposition ou le maximum dintensité augmente indépendamment, la hauteur et la pente de Vonde-b augmentent 
tous deux et tendent a la saturation. Le temps de latence ne montre que peu vu pas de changement. La relation de BUNSEN-ROSCOE 
entre le flux maximum et le temps d’exposition semble se maintenir pour un observateur. 

Quand le temps d’établissement est changé, le temps d’exposition et le maximum d’intensité restant constants, aucun change- 
ment rapparait dans VERG bien que Vénergie totale du stimulus ait été doublée. L’analyse de VERG cherchant a vérifier la 


relation BUNSEN-RoscoE est de ce fail considérablement simplifiée : 


une simple intégration. 


le mécanisme des bdtonnets ne peut étre considéré comme 


Des stimuli de breve durée donnent la réponse de Vonde-x (caractéristique du mécanisme des cénes) qui se manifeste de trois 


facgons différentes a des fréquences sensiblement égales. 


L’analyse des résultats obtenus pour des stimuli standards sur une période de plusieurs mois montre une variation saison- 


niéere typique et appréciable dans la réponse de Vonde-b. 


ZUSAMMENFASSUNG. 


Am Menschenauge wird die Abhdngigkeil des ERG von der Lichverteilung tn einem weissen Reizlicht wah- 


rend der Dunkeladaption bei zwei Beobachtern gemessen. Der Reiz wird nach drei Verdnderlichen aufgegliedert : Anstiegszeit, 
Gesamtdauer und Gipfelwert. Die erhaltenen Kurven werden nach drei charakteristischen Werten analysiert : Hohe, Steigung 


und Latenzzeit der b-Welle. 


Wenn man unabhdngig voneinander die Gesamtdauer oder den Gipfelwert erhéht, so nehmen auch Héhe und Steigung der 
Welle zu und ndhern sich einem festen Sdttigungswert. Die Latenzzeit zeigt hingegen keine oder nur eine unbedeutende Ande- 
rung. Das BUNSEN-RoscoE sche Gesetz zwischen Gipfelwert und Gesamtdauer scheint fiir etnen der Beobachter zu gelten. 

Bei einer Anderung der Anschwellzeit des Reizes unter Beibehalten der Gesamtdauer und des Gipfelwertes zeigte sich keine 
Anderung des ERG, obwohl die Gesamtenergie des Reizes auf das Doppelte erhoht wurde. Die Analyse des ERG im Hinblick 
auf eine Bestdtigung des BUNSEN-RoscoE schen Gesetzes ist deswegen ein allzu einfaches Vorgehen. Der Stébchenmechanismus 


kann nicht als ein einfacher integrierender Vorgang angesehen werden. : ; ae ; . : 4 
Kurzzeitige Reize ergeben die x- Welle (charakteristisch fiir den Zapfenmechanismus), die in dret Arlen ndherungsweise gleicher 


Hadufigkeit auftritt. 


Die Analyse tiber Monate sich erstreckender Messungen mit Normatreizen zeigt in der Hohe der b-Welle eine typische jah- 


reszeitliche Schwankung von betrdchtlicher Grosse. 


Introduction. — From time to time parallels have 
been drawn between some characteristics of human 
vision and the electrical response to light for the eyes of 
the lower animals, with a word of caution with regard 
to the difference in spevies. Now, with the advent of 
satisfactory contact-lens techniques in the recording 
of the human electroretinogram (ERG), experiments 
are being conducted in some laboratories in both the 
objective and subjective fields. 

Attemps to correlate the electrical and psychophy- 
sical responses of the eye take on, therefore, a greater 
significance. However, the removal of the earlier 
obstacle of a species difference only underlines even 
more strongly the necessity of proceeding carefully 
in the drawing of parallels. 


The present work is part of an investigation into 
the effect on the ERG of the time distribution of a 
light stimulus. The growth and decay of the stimulus 
employed is shown in figure 1. The rate of rise (over 
section A B) and the time of exposure (section A D) 
are independently variable. In an initial study repor- 
ted by one of us [1] the duration of the steady level 
BC was made effectively infinite (i. e. greater than 
0. 6 sec. for recording purposes), and only the time of 
rise was varied. This type of stimulus was called a 


* The research reported in this document has been made 
possible through the support and sponsorship of the Air 
Research and Development Command, United States Air 
Force, through its European Office, under contract No. AF61 
(514)-634 C with the Istituto Nazionale di Ottica, Arcetri- 
Firenze (Italy). 
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Fic. 1. — Generalised stimulus form. 


«time-gradient of illumination ». Results were obtained 
for both light- and dark-adaptation. It was found 
that both the shape and size of the ERG were inde- 
pendent of the time of rise if this was less than 70 msec. 
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Ig, 2. — Falling plate on objective aperture. 
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For higher values the amplitude of the b-wave remai- 
ned constant and its slope decreased while the latency 
time increased. 

The present paper reports experiments in dark 
adaptation using a wider variety of stimulus forms in 
which time-gradients of illumination occur as a limi- 
ting case. 

Some interest for this paper rests on the BUNSEN- 
Roscoe law. Granit [2] has reviewed the connec- 
tion of earlier electroretinographic work with this 
subject. 

In a more recent study Jounson and Barrett [3 | 
have investigated the effect of duration and intensity 
of a ight stimulus on the b-wave of the human ERG. 
It was found that in order to produce a b-wave of 
some fixed amplitude, the intensity and exposure 
time of the stimulus were interchangeable for expo- 
sure times less than 100 msec. For times greater than 
this only the intensity played a part. 


I. Apparatus. — The apparatus has been described 
in an earlier paper [1] but there have been some 
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Fic. 3, — Variation of height (/y,), slope (s},) and 
latency time (t,) of the b-wave with exposure time 
({;) for a constant peak intensity. 
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changes in conditions. A summary is given below 
with an amplification of some points relevant to the 
present paper. 


The Stimulation System. A white light source is 
imaged by an objective on to a 3 mm diameter circu- 
lar aperture. This aperture is viewed through an eye- 
lens, which projects it to infinity and defines the field 
of view, which is 3.49, approximately. The eye-lens 
images the objective at the eye pupil, but since the 
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objective aperture is used to control the stimulus 
intensity, the size of this image is variable. However, 
even for the largest aperture this is not more than 
2 mm. 

A slot P (fig. 2) in a plate falling under gravity cros- 
ses the light beam close to the objective and passes 
the light stimulus. The objective is stopped down to 
the shape of the trapezium Q (for a linear rate of 
rise of intensity, Q should be provided with sides 
which are suitable parabolic arcs). 


M. (d.a.) 
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Variation of the ERG with exposure time (f,) for a constant peak intensity. 
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The time of rise of the stimulus depends on the in intensity occur simultaneously on all parts of the 
drop-height of P and the height J of Q, while the drop- _ retinal image, a requirement fundamental to the pre- 
height of P in combination with the length L of P sent study. 
determines the time of exposure. The stimulus inten- 
sity is controlled by the area of the trapezium. The Recording System. — Electrical potentials are 

This optical system ensures that the same changes picked up between a silver electrode mounted on the 
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Ilia. 5. — Variation of height (hy), slope (s},) and latency time ({,) of the b-wave with exposure 
time (f,) for a constant peak intensity, using three times of rise (tp). 
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corneal bulge of a contact lens and an indifferent 
silver electrode on the forehead. The signals are pas- 
sed through a capacitance-coupled amplifier to an 
oscilloscope. At the amplifications which are used 
the noise level is about 10 uV. The sweep of the oscil- 
loscope is provided by means of a circuit activated 
by the falling plate P. The responses are recorded 
photographically simultaneously with the form of 
the light stimulus. 


II. Method. — After inserting the contact lens the 
observer dark adapts for some 20 minutes or so. The 
head is steadied by the dental impression method 
and the test field fixated foveally. In the absence of the 
stimulus, fixation is maintained by a red spot of low 
luminance. Between the presentation of each stimulus 
and the next, a period of dark adaptation is allowed. 


This interval varies from 3 mn to 5 mn depending 


on the duration of the stimulus employed. Each sti- 
mulus is presented at least 10 times and in each 3es- 
sion the responses to a standard stimulus also are 
recorded to keep a check on the constancy of condi- 
tions and day-to-day variations. 

On one or two occasions a general drift was found 
in the magnitude of the response to the same stimu- 
lus throughout a session. This could not always be 
attributed to incomplete dark-adaptation, but such 
records are excluded from the present analysis. 

Two observers took part, J. M. and L. P. Over a 
period of some months a large general drift in the 
response to the standard stimulus was noticed (see 
Appendix). 


III. Symbols. 


MCUSTIMUlUS, — Lime Ol LISCIM MISCO... 246 ese: ie 
Time of exposure in msec ...... t 


The b-Wave Response. — Maximum height in pV. hy, 
Rate of riseinuV/msec.  s; 
Latency time inmsec.. f¢ 


IV. Experiment 1. — The maximum level of the 
stimulus is kept constant at 100 lux at the pupil 
plane. 

This statement is rather ambiguous when it refers 
to short duration stimuli. It is understood here as a 
specification of the luminance of the field if the hori- 
zontal section BC (fig. 1) were continued for an inde- 
finite period. 

Keeping f) constant at 10 msec., records were 
obtained for t,, varying from 20 msec to infinity (i. e. 
greater than 600 msec). In figure 3 are shown the 
variations of h,, s,, and t, with changes in ¢,. It is 


Fic. 6. — Variation of the ERG with the time 
of rise (tj), keeping peak intensity constant, and for 
time of exposure ((,) effectively infinite. 
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seen that f,, increases with ¢, rapidly at first and 
then more slowly, tending towards saturation. A 
similar change occurs for s,, but ¢, remains unaffected. 

The response curves for J. M. are shown in figure 4, 
each being the mean of 10 or more records as are all 
the results reported below. The average deviation of 
each set of records for any particular session is 
about 20 wV over the range of the trace which is 
considered in the present analysis. Below each trace 
the form of the light stimulus is indicated. Here we 
see that hy, increases with ¢, for stimuli with ¢, < ¢, 
(t, being constant at 50 msec). But when ¢, > ¢, then 
hy, no longer depends on t.,.. 

To verify whether or not the time of rise played 
any part here, the experiment was repeated for 
t= 2 msec and for 7, =30) msec. The results for 
L. P.are shown in figure 5. The same has been verified 
for J. M. The values for ¢) = 10 msec are those found 
in figure 3. No points are possible for ¢, = 30 msec 
when ¢,, is less than 60 msec since the stimulus form 
is symmetrical. 

It is seen that within experimental error the three 
sets of observations coincide. The same agreement is 
found for J. M. 

It is useful to recall at this point some results 
obtained previously [1] for another observer S. G. 
in which time-gradients of illumination were employed 
using a 1° field but the same peak intensity. These are 
summarized in figure 6. 
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70 <t, < 140 msec the value of sj, fell, i, remaining 
constant, and for t) < 140 msec a tendency for both fy, 
and s), to decrease was shown. There was also a sug- 
gestion that ¢, rises slightly for f > 100 msec. 


V. Experiment 2. — It was suggested that an 
increase in the peak intensity of the stimulus might 
compensate for a reduction in exposure time if the 
fractional changes in these two variables were the 
same. With this question in mind one of the stimuli 
used in Expt. 1° (t) = 410 msec, ¢,, = 20. msec) was 
chosen and records were obtained for relative peak 
intensities of 1, 2, 3, 3.5 and 5. The average traces 
for J. M. are shown in figure 7 in which the stimulus 
forms are also indicated. In this experiment variations 
in the rate of rise of the stimulus must also take place, 
but it has been shown (see fig. 5) that changes in the 
time of rise from 30 to 2 msec for a given peak intensity 
level produce no change in the ERG. It seems 
improbable that a serious disturbance is introduced. 

In figures 8 and 9 for subjects L. P. and J. M. res- 
pectively the traces have been analysed as before 
into variations in h,, s, and ¢,. For both subjects it is 
found that an augmentation in peak intensity produ- 
ced first an increase in fj, and later no appreciable 
change. 

The same kind of change for s, is recorded for 
subject L. P. but not for J. M. for whom s, shows 
no tendency to saturate. Latency time ¢, shows a 


It was found that in each case for tj << 70 msec the small diminution for L. P. but remains constant 
b-wave was unaffected while in the interval for J. M. 
pv 
300 
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‘Ia. 7, — Variation of the ERG with increasing peak intensity for constant exposure time ((,,). 
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Fic. 8 and Fria. 9. — Variation of height (hj), slope (s}) and latency time (/,) 
of b-wave with increasing peak intensity. 


In order to test the possibility of a reciprocal rela- 
tion between peak intensity and.exposure time the 
results of experiments 1 and 2 are compared. The 
method of comparison may be best explained by 
reference to figures 10 and 11. 

The ratios L/E, are plotted as abscissae where FE 
is the total energy contained in the stimulus conside- 
red while F, refers to the standard stimulus for 
which ¢t,, = 20 msec and where the relative peak 
intensity is 1. As ordinates the ratios h/h» are plot- 
ted where h and /y are the heights of the b-waves cor- 


eats 


3} ho 


responding to £ and £, respectively. The broken 
curves derive from Experiment I in which F is raised 
by increasing ¢,, and the full line curves are obtained 
from Experiment 2 in which F£ is raised by increasing 
the peak intensity of the stimulus. 

Mention should be made here of the derivation of 
the ratio h/hy. It was found early on in the experi- 
ments that variations in the response to a standard 
stimulus occurred from session to session. If these 
variations are part of some general change in the 
observer’s condition (see Appendix) then their 


Varying. te 
Varying | (peak intensity) 


h 


Subject J.E 


Fic. 10 and Fic. 11. — Relative changes in height of b-wave (/),) with increasing exposure time (f,,) 
at constant peak intensity and with increasing intensity (J) at constant exposure. 
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effect in the present analysis could be reduced by 
taking values of 2 and Ay which were obtained in the 
same session. 

For observer L. P. it is seen that in order to pro- 
duce a given value of 2, it is only necessary that the 
energy contained in the stimulus should be constant, 
at least for the range of variations that have been 
considered. The same cannot be said for observer 
J. M. for it is clear that changes in intensity at cons- 
tant exposure time are more effective than increases 
in exposure times at constant intensity. 

The divergence between the two observers is con- 
firmed also in the following way : two points A and B 
are considered on the curve in figure 5 showing the 
variation of Ay with ¢, for constant peak intensity. 

For A, ¢, = 60 msec and for B, ¢, = 20 msee. On 
reducing the intensity for A to a third, the new b- 
wave was recorded and found on subject L. P. to be 
equal to that for B but not so for J. M. Reducing 
intensity also involved a change in the rete of rise of 
the stimulus, but it has already been shown that this 
plays no part. 


13. — Modes of the x-wave response. 


VI. Experiment 3.— For the stimuli which are 
used in previous experiments no definite evidence 
of cone activity is recorded, i. e. both the a-wave 
and the x-wave are either absent or barely distin- 
guishable above the circuit « noise » level. 

According to Wrretsa [4] the cone contribution is 
brought into evidence a the time of exposure 
is less than 20 msec. 

Using the stimulus ¢, © 
find that the x-wave emerges clearly. 

Of the responses examined (totalling 50 for the 
above stimulus) there appear to be three forms which 
occur with equal frequency (fig. 12). In this figure 
each response is spht mto two waves more or less 
separated one from the other but which, from their 
latency times, may be identified as the x-wave and 
the b-wave. 

Two more stimuli for which ¢, is 20 msee but for 
which tf changed from 2 msee to {0 msee were exa- 
mined. In each of these cases the x-waves found were 
usually very small, in fact hardly perceptible above 
the background «noise» of the cireuit (fig. £8). 


= 4 msec, t, = 2 msec, we 
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However, the frequency of occurrence of the x-wave 
is greater for the shorter time of rise (see table 1). 


Pare f 
Stimulu %-Wave: Y, Re ponse| Total No® 
of 
tyngee t, Insc Absent Present | Response 
| ¢ 

A 2 4 ) 100 50 
B Zz 20 33 66 G2 
C 10 20 65 15 104 


VII. Discussion. In the first two experiments 
the findings of Jounson and Barriers [ 3} mention- 
ed above have been partially confirmed for one sub- 
ject L. P. but not for the other J. M., and for L. P. a 
BuUNSEN-RoscoE relation for stimulus peak intensity 
and exposure time seems to hold. However, the cri- 
tical value for exposure time (i. e. the value of ¢ 


f 
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above which no change in hy is recorded) is found to 
be of the order of 50 ~ 60 msec as compared with the 
value 100msec obtained by Jounson and Barrett ; 
but since our viewing conditions are different, per- 
haps this is not surprising. But on Jooking again at 
Experiment 1 it will be clear that it is not possible 
to speak of the energy of a stimulus independently 
of its form, for it will be remembered that when the 
effect of varying the time of rise of a stimulus is tes- 
ted, keeping the time of EX POBUSE and peak intensity 
constant, it is found to produce no effect either in h,, 
s, or t, although the total energy contained in the 
stimulus might be halved and although the energy 
might be delivered with the latency period. This 
indicates that the carrying the BuNseN- 
toscor relation, from psychophysical threshold 
studies to the electrical response of the retina, is an 
oversimplification. This is illustrated perhaps more 


clearly in figure 14. 
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In this figure are summarized the results from two 
observers : L. P. and J. M. and another observer 
reported in an earlier paper S. G. [1]. In the columns 
are shown equivalent stimuli, i. e. stimuli producing 
the same ERG (no change in h,, s), or t,). Columns a, 
b, ec, d, refer to experiments 1 and 2 in which row 
A applies to both observer L. P. and J. M. in expt. 
1, while row B and C apply only to L. P. in expt. 
2. Glancing down any one column brings out clearly 
that while it is possible to find two stimuli of the 
same energy which produce the same response, it 1s 
also possible to find others of different energies which 
do the same. 

Alternatively looking from column to column we 
find stimuli of the same energy which produce diffe- 
rent responses. In row B the effect of increasing 
intensity is shown to produce an increase in //), up 
to t,, = 30 msec, beyond this hj), remains constant but 
additional energy seems to be used in producing a 
shorter latency time. 

Row A, column e, refers to the earlier work on 
subject S. G. on time gradients of illumination. Here 
it was found that fj, is independent of the stimulus 
form but that when the time of rise exceeds the 
latency time the value of s,, falls off. 

In the first two spaces and column a of row A are 
shown stimuli which produce a mixed rod-cone res- 
ponse. The shorter the exposure time, the greater 
the frequency of occurrence of the x-wave. 

Our results seem in agreement with the fact shown 
by all the work on sensation from the days of KONIG 
up to the present time, and also the objective rese- 
arches based on the electroretionographic response, 
that the dark adapted eye is highly sensitive but 
slow, like a ballistically recording galvanometer, and 
it has been assumed to integrate the total quantity 
of energy reaching it [5]. 

However, this cannot explain the apparent loss of 
energy that we have shown for rapidly rising stimuli, 
and as has been mentioned above, the total energy 
contained in a stimulus cannot be considered inde- 
pendently of the form. 

We intend to continue our study in the light adap- 
tation state. In this case, according to Wricut [6] 
and Granit [7], the response should consist of very 
complicated excitation inhibition patterns giving 
the higher centres a cue for discrimination. 


Appendix. — It was noted early on in the work 
that variations in the height of the response toa 
given stimulus occurred from day to day and, as 
mentioned above, this could be allowed for, where 
necessary, by taking ratios. After a period of some 
months, however, it was found that a general drift 
in the size of the response was taking place above the 
day to day fluctuations. 

In the earlier work [1] only the stimuli, termed 
time-gradients of illumination, were used and_ since 
the oscilloscope employed had been changed twice 
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Fra. 15. — Seasonal variations in the height of 
the b-wave (/),) for a standard stimulus. 


during that period, the change in response was attri- 
buted to changes in circuit conditions, and the varia- 
tion in the standard stimulus (¢, = 10 msec) was 
used as a correction factor for the calibration. These 
records refer to the period August 1955 to February 
1956. For the later work, reported in the present 
paper, the standard stimulus used was a symmetrical 
flash with ¢,, = 20 msee and ¢, = 10 msee. Both this 
and the earlier used standard stimulus gave rise to 
characteristic rod responses. 

The variations in the height of the b-wave of this 
response are shown in figure 15. Each point is the 
mean value of hj, for the standard stimulus employed 
in any one session. 

The variations are seen to be considerable over a 
period of months. Allowing for observer differences 
and assuming continuity over the whole period, 
there is evidence for a variation which is typically 
seasonal : i, e. maximum and minimum activity 
occurring in summer and winter respectively. 
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An uncertainty relation in interferometry 


B.S. Tuornton 


Division of Metrology, National Standards Laboratory, Commonwealth Scientific 
and Industrial Research Organization, Sydney, Australia. 


In multiple-bear interference microscopy the measu- 
rements made of surface details can be of almost mole- 
cular magnitude, and it might be expected that the 
HEISENBERG uncertainty principle for photons will 
set some limit to the accuracy of measurement of 
surface details revealed. 

Starting with the relation Ak, Ax~| for the sprea- 
ding of wave groups, where k,, is the z-component of 
the wave number vector, InceLsram [1] observed 
that Z/N was the sine of the total variation of the 
wave number vector in the Az region of the wave 
field of the interferometer, N being the order of inter- 
ference, Z the mean direction of propagation of the 
light and % the wavelength used. The spreading rela- 
tion given above then gives 


(1) Ag Az ~ Ni? 


InGeLstam gave V2? only as an order of magnitude, 
and for the exact relation writes 


(2) Az Az > = N2? 


where «/8 is a constant (unspecified) which, amongst 
other things, depends on the order of interference. 

An independent approach to the uncertainty pro- 
duct (2) can be obtained as follows. The fundamental 
equation of interferometry is 


(3) Ni = 2 nt cos 8 


2nt cos@ being the optical path difference between 
successive interfering waves, and r being the refrac- 
tive index. Since light is no longer strictly monochro- 
matic if it is not emitted over an infinite time interval, 
and since light sources in practice emit truly coherent 
wave trains only for periods of the order of 10-7 or 
10-* seconds, it is necessary to consider a wave group 
which, as can be shown by Fourier’s theorem, is 
equivalent in all its effects to the superposition of 
wave trains of different wavelength ranging from 
1 > to 4,, + ou , where 2,, is the mean wave- 
length and Aa the range of wavelengths present. In 
optics the wave groups are of quite large dimensions 
‘compared with the mean wavelength 4,,. From Equa- 
tion (3), 


(4) s=>= 


(since ¢ is measured in the z direction). The momentum 
p of a photon has the same direction as the propagation 
vector k and has a magnitude p = h/2 so that conside- 
ring magnitude only, 


(5) Ap = a 
P A 


From Equations (4) and (5) 
h hz 


A 


and the Hrtsenperc Uncertainty Relation states 
that : 

(7) Ap, \ x > h/4n 
so that from equations (6), (7) and (3) 


N22 


8 Lh 
(2) “8 an cos 0 


which corresponds to the relation (2) with «/8 equal 


to : 
8xzncos 6 


That there is an uncertainty in z can also be seen 
by considering the electromagnetic field Bey eee) 
incident on an interferometer with the appropriate 
boundary conditions to determine the fringe pattern. 

The intensity distribution of the emergent light is 
given by 

7 v* %* 

(9) a EE* + HH 

8x 


where / and H are complex amplitudes determined 
from the boundary conditions for the system, the 
starred quantities denoting the complex conjugates. 

The number of photons corresponding to J is given 
by 

EE* + HH* — 
where A(z) = 1/k, 
and for a region large compared with 4, P(x) is ana- 
logous to ¥ ¥* in the quantum mechanics of electrons 
and other particles. 

The HetsenperG Uncertainty Relation for the elee- 
tromagnetic field [2] states that : 

he 


hia 
2 4l 8p 


(10) 
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where Sv = (81)? is a small element of volume whose 
value depends on the method of measurement used, 
and c is the velocity of light. 

The uncertainty in J and P(z) arising from (10) 
gives rise to an uncertainty in z through the relation 


(11) AP(z) AD W14 


between the phase, ®, of the wave and the number of 
photons [3] P(x) where, in the present case 


AD = 2m Aza. 


In practice the value for Az is governed by the sharp- 
ness of the fringes since a fringe can be located to 
within only a certain fraction of its half width. 

Essentially, equation (8) means that when a very 
great accuracy in depth measurement is required an 
unlimited accuracy in surface resolution cannot be 
attained. 

ToLansky and Emara [4], using low order mul- 
tiple-beam fringes, have been able to use a multiple- 
beam interference microscope with high magnifi- 
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cation (x 1500) to resolve volumes of the order of 
5 x 10-4 ce which, according to the uncertainty 
relation (8), is almost the limit obtainable, since the 
minimum resolvable distance in practical microscopy 
is known to be 4/2, so: 


Nw 


Az Aa Ay = ———————— 
ses Smn cos 0 


r 
2 


i. e. the minimum resolvable volume is © 3.5 « 10-1 ce 
for V =1,n061, 6S 0,X%5.5 < 10" om: 
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Quelques propriétés optiques des milieux stratifiés symétriques 


F. ABELES 
Institut d’Optique, Paris 


Nous indiquerons, dans ce qui suit, quelques pro- 
priétés des milieux stratifiés symétriques, qui n’ont 
pas encore été signalées, tout au moins a notre con- 
naissance. Un milieu stratifié est limité par deux plans 
paralléles et son indice de réfraction n —jk = Ve 
(réel ou complexe) et sa perméabilité magnétique pu 
sont constants dans des plans paralleles a ces deux 
plans et ne varient que dans la direction OZ perpen- 
diculaire 4 ces plans. Le milieu stratifié sera dit symé- 
trique lorsque e(d —z) = ¢(z), w(d — z) = u(z), ad 
étant son épaisseur totale, étant entendu qu il 
s’étend entre les plans z = 0 et z = d. 

Nous supposerons d’abord que le milieu stratifié 
symétrique est absorbant (k= 0), et qu’il est limité 
par deux milieux non absorbants homogeénes d’indices 
de réfraction n, et n,. Nous désignerons par r, r’ et t 
les coefficients (complexes) de réflexion pour une 
onde incidente du coté du milieu d’indice ny ou n, 
et pour une onde transmise respectivement (1). Les 
facteurs de réflexion et de transmission seront dési- 
gnés par R, R’, T, tandis que les déphasages par 
réflexion et par transmission seront 5,, 8,, et 3,. Nous 
aurons donc r = R"? ei, etc. Expérimentalement, 
on mesure, pour une incidence o, et une vibration rec- 
tiligne perpendiculaire (1) ou paralléle (||) au plan 
dimerdence, ART, «. = 8.-— 2)V yy, 1) =o, 
é, = 8, + You, ot Yo = my COS g et y =(2 @ ji. 


(1) Rappelons que ces « coefficients » sont des rapports d’am- 
plitudes (complexes) et contiennent donc aussi les déphasages. 


1. Dans le cas d’un milieu stratifié symétrique 
absorbant, ces six quantités mesurables ne sont pas 
indépendantes. Il y a entre elles la relation suivante : 


Pees 
(—re—m—2|0—m Vee cos «7 + 
0 


+T 
+ 2V RR’ Tcos (e+ © p— 2 3) 0, 


ou Ry représente, pour la méme incidence @, et la 
méme polarisation, le facteur de réflexion lorsque les 
milieux my) et n, sont au contact (c’est-a-dire lorsque 
d =0, d’ot Vindice 0). 

La démonstration de cette relation sera donnée plus 
tard, dans un mémoire plus développé. Disons ici 
qwil était a prévoir qu'il doit y avoir une relation 
entre les six quantités mesurables, car la matrice 
caractéristique du milieu stratifié symétrique ne con- 
tient que cing paramétres indépendants. 


2. Pour une incidence @ donnée, les six quantités 
mesurables relatives A une vibration | ou |] per- 
mettent de déterminer l’épaisseur du milieu stratifié 
symétrique, sans que les lois de variation de e(z) et de 
u(Z) nous soient connues. Pour cela, il suffit de résoudre 
Péquation suivante : 


asin2 Yon + beos2 Yyn +c=0 
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ou 


7 R 
a =—\/RR cos(e,-+ ¢.7)+ vee cos ¢,+ T cos 2 «, 
9 


age: a R 
b=V RR’ sin (e, + a sin ¢,—T sin2 ¢, 
0 


Remarquons que hous avons pose 
Vey gl Gra kf) 


dans le cas de l’incidence normale, les n,; étant 
remplacés, suivant le procédé habituel, soit par 
Yi; =n, cosg, soit par Z; = n/cos 9;, selon 
quil s’agit dune vibration | ou || respectivement 
(n, SIN 9, = Ny sin »® détermine l’angle de réfraction 
dans le milieu n,). 

Nous laisserons de cdté ici une discussion plus 
détaillée de cette relation. 


3. Lorsque le milieu stratifié symétrique se réduit 
a une seule couche mince absorbante, il y a, de plus, 
entre les six quantités mesurables relatives a la vibra- 
tion | et les mémes quantités relatives A la vibra- 
tion ||, Pangle d’incidence ¢) étant le méme dans les 
deux cas, la relation suivante : 


Toy legal ad To\| R 


x x = 
Roy e Ral T\ 


1B 
a To fps 
ty [Fonte | / i Rpc0s 61, 24 — ey re) 


; 1 
ee cheat — Seu > Fx Ey) Ene 


ou 7, est l’analogue, pour londe transmise, de Roy 
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défini plus haut. La convention de signe relative 
a VR, est toujours celle quia été faite plus haut. 


4. Dans le cas d’un milieu stratifié symétrique non 
absorbant (k = 0), il n’y a qu’un facteur de réflexion 
et un déphasage indépendants. Entre les six quanti- 
tés R, R', T, 8,, 8,7, et 3, il y a les quatre relations 
suivantes : 


RSRS1—7T; 3-3 293, So, 


cos 8, + ce cos (8, — 8,) = 0, 
0 


la derniere relation étant indiquée ici pour la premiére 
fois, vraisemblablement. Elle permet, connaissant, 
par exemple, R, «, et ¢, (quantités mesurables) de 
déterminer lépaisseur du milieu stratifié symétrique 
non absorbant, sans connaitre les lois de variation 
de «(z) et de u(z), a partir de la relation 


R 
cos e, + R, °% (e,— €,) 
0 


tg Yon = — 


. ie 
—— Sine R, sin (¢,—€,) 


Ceci montre aussi que pour un angle d’incidence 
donné, la fraction qui se trouve a droite du signe égal, 
est la méme pour une vibration | et pour une vibra- 
tion || au plan d’incidence. 

Nous ne discuterons pas ici les applications pra- 
tiques de ces résultats. Remarquons simplement que 
les milieux stratifiés symétriques sont tres souvent 
utilisés dans la pratique (par ex. les systemes de 
lames minces alternées comportant un nombre impair 
de couches de haut et bas indice et les filtres interfé- 
rentiels du type Fasry-PERor simples ou com- 
plexes). 


Manuscrit regu le 16 octobre 1956. 
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The interference systems of crossed diffraction gra- 
tings (Theory of Moiré fringes). J. Guirp. The 
Clarendon Press, Oxford, 1956, pp. 152+ vi. 
25 s. net. 


The moiré fringes resulting from crossed diffraction 
gratings, noted by Ray.ereH in 1874, have recently 
been used by Merron as a test for the ruling errors in 
his gratings. Arising from Merron’s elegant work, the 
National Physical Laboratory has for the past seven 
years undertaken a programme of development of 
grating production and testing, and an important 
sideline of this work has been the application of the 
moiré fringes themselves to automatic metrology. The 
fringes resulting from two replicas of the same gra- 
ting placed almost in contact can be made to indi- 
cate relative motion between the replicas ; and the 


possibility of counting the fringes photoelectrically 
has opened up the prospect of automatic machine 
tool control, by replacing the normal machine scales 
by grating replicas. 

Mr Guitp has been responsible for much of the 
work at the N. P. L. on the testing of gratings and it is 
therefore of considerable general value that he should 
have produced the present book, which describes his 
development of the theory of moiré fringes and illus- 
trates the influence of grating defects. He has written 
‘“‘ with the primary object of helping those engaged in 
harnessing moiré fringes in the service of technology 
and industry ” and the book is concentrated on the 
fundamentals of the moiré patterns themselves. The 
effect of ruling errors on the shape of the fringes, and 
the influence of groove form on intensity distribution 
are the author’s main concern in the present mono- 
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graph. He has reserved technical matters, such as the 
best ways of observing the fringes and of mounting 
the gratings, for a shorter paper which is in prepara- 
tion. 

This is not an easy book to review, in that one needs 
to acquire a working familiarity with the nomencla- 
ture which the author has had to develop for this new 
field, and the detail, although straightforxward, is 
necessarily heavy. A superficial comment is that the 
reader would probably be grateful for more ‘‘ signposts”. 

The author, for example, limits his treatment of 
the “‘ crossing ”’ of gratings for most of the book to 
instances where the grating surfaces are paralled, and 
the reader inexperienced in this field may ask himself 
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why this limitation is imposed, when in practice the 
surfaces of any two gratings are unlikely to be exactly 
paralled. The author adequately answers this question 
at the end of the book (and indeed the answers to 
many questions concerning fringes can be found in 
the book) but an anticipatory reference would have 
been helpful at the beginning of the analysis. 

Although there are few diagrams in the text, there 
are some good photographs at the end of the book 
showing the characteristics of the principal kinds of 
fringe defect. Mr Guitp has given an impressive 
account of this aspect of the work behind the applica- 
tion of moiré fringes. 


R. V. Jones. 


INFORMATION 


Colloque de spectroscopie interférentielle. — Un Col- 
loque international, organisé par le Centre National 
de la Recherche Scientifique aura lieu aux Laboratoires 
de Bellevue, pres de Paris, du 9 au 13 Septembre 
1957 sur le sujet suivant : Progrés récents et méthodes 
nouvelles en spectroscopie interférentielle. 

Suivant la formule adoptée par le C. N. R. S. pour 
les Colloques Internationaux, le sujet du colloque est 
nettement limité et la participation au colloque est 
réservée & ceux qui ont contribué au développement 
du sujet par leurs recherches personnelles ; cependant 
les physiciens qui sont spécifiquement intéressés par 
le sujet peuvent assister au colloque a titre d’auditeur. 

Plusieurs physiciens d’Europe et d’Amérique ont 
déja promis d’assister au colloque et d’y présenter 
des communications. Entre autres : P. FELLGeErt, 
D.A. Jackson, J.H. Jarret, F. A. Jenkins, H. Kunn, 
J. E. Mack, K. W. Meissner, J. Ranp Mac Natty, 
A. STEUDEL, J. Stronec, J. TERRIEN. 

Les communications prévues sont, provisoirement, 
groupées sous les rubriques suivantes : 


I. Principes généraux des méthodes nouvelles de 
spectroscopie interférentielle (par exemple : Fasry- 
Perot photoélectrique, méthodes par modulation et 
transformation de Fourier, etc.). 


II. Problemes de principe spécifiques au Fasry- 


Perot (par exemple : extensions aux domaines spec- 
traux étendus, aux problemes astronomiques, FaBRy- 
Perot sphérique, etc.). 


III. Procédés d’exploration des spectres avec le 
Fasry-PEror (variations de pression, déplacements 
mécaniques, utilisation de la magnétostriction, etc.). 


IV. Qualités des lames et couches réfléchissantes 
pour le Faspry-Prerort ; détails pratiques. 


VY. Descriptions et performances d’équipements 
Fapry-PeEror. 


VI. Problémes métrologiques. 


VII. Mise en pratique des méthodes interféren- 
tielles autres que Fasry-Peror. Instruments et résul- 
tats. 


Cependant toute contribution qui n’entre pas dans 
les paragraphes précédents pourra étre acceptée a 
condition qu’elle soit nouvelle et intéresse les méthodes 
de spectroscopie interférentielle (les propriétés et les 
modes classiques dutilisation des réseaux sont exclus 
du champ du colloque). 

Des renseignements complémentaires peuvent étre 
obtenus en écrivant 4 Vadresse suivante: Professeur 
P. Jacquinor (Colloque S. I.), Laboratoires du 
C. N. R.S., Bellevue (S.-et-O.). 
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